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ABSTRACT 

We treat of the high-energy astrophysics of the inner ~200 pc of the Galaxy. Our 
modelhng of this region shows that the supernovae exploding here every few thousand 
years inject enough power to i) sustain the steady-state, in situ population of cosmic 
rays (CRs) required to generate the region's non-thermal radio and TeV 7-ray emis- 
sion; ii) drive a powerful wind that advects non-thermal particles out of the inner GC; 
iii) supply the low-energy CRs whose Coulombic collisions sustain the temperature and 
ionization rate of the anomalously warm, envelope H2 detected throughout the Cen- 
tral Molecular Zone; iv) accelerate the primary electrons which provide the extended, 
non-thermal radio emission seen over ~150 pc scales above and below the plane (the 
Galactic centre lobe); and v) accelerate the primary protons and heavier ions which, 
advected to very large scales (up to ~10 kpc), generate the recently-identified WMAP 
haze and corresponding Fermi haze/bubbles. Our modelling bounds the average mag- 
netic field amplitude in the inner few degrees of the Galaxy to the range 60 < B / fiG 
< 400 (at 2(7 confidence) and shows that even TeV CRs likely do not have time to 
penetrate into the cores of the region's dense molecular clouds before the wind removes 
them from the region. This latter finding apparently disfavours scenarios in which CRs 
- in this star-burst-like environment - act to substantially modify the conditions of 
star-formation. We speculate that the wind we identify plays a crucial role in advect- 
ing low-energy positrons from the Galactic nucleus into the bulge, thereby explaining 
the extended morphology of the 511 keV line emission. We present extensive appen- 
dices reviewing the environmental conditions in the GC, deriving the star- formation 
and supernova rates there, and setting out the extensive prior evidence that exists 
supporting the notion of a fast outflow from the region. 

Key words: cosmic rays - galaxies: star formation - Galaxy: center - ISM: jets and 
outflows - ISM: supernova remnants - radiation mechanisms: non-thermal 



1 INTRODUCTION 

In recent years a picture has emerged that the inner regions 
of star-forming galaxies should i) be important sources of 
7-rays in the universe (Pavlidou & Fields 2002 Thompson 



et al.|2007| [Dogiel fc Breitschwerdt|2009[); ii) drive powerful 
galactic winds ( Dogiel fc Breitschwerdt|2009 1 and iii) there- 
fore, be important shapers of the inter-galactic medium, par 
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ticularly its metallicity ( Strickland fc Heckman||2009 K The 
Galactic centre (GC) presents our closest view of a Galactic 
nucleus. The physical conditions pertaining in the GC put it 
on an intermediate footing between the relative quiescence 
of the Galactic disk interstellar medium (ISM) and the envi- 
rons at the centre of luminous starbursts. Still, when viewed 
from the local perspective, many of the GC ISM parame- 
ters are extreme. The GC proffers, then, perhaps our best 
hope for a detailed understanding of the processes shaping 
the evolution of the nuclei of spiral galaxies and, in general, 
regions energised by galaxy mergers or interactions ( )Melia| 
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2007). Such processes include star formation in giant molec- 



ular cloud complexes, concomitant supernova activity, the 
resultant generation of nuclear outflows or winds, and the 
acceleration of non-thermal particle distributions. 

This paper is the third in a series related to the high- 
energy astrophysics of the Galactic centre (GC). In the first 
paper of this series ( Crocker fc Aharonian^, 2010, Paper I) 
some of us showed that it is likely that star-formation and 
consequent supernova processes in the inner ~200 pc of the 
Galaxy have, over timescales approaching the age of the 
Galaxy, injected cosmic ray (CR) ions to large heights above 
and below the Galactic plane creating thereby the very ex- 
tended regions of diffuse, non-thermal emission recently un- 
covered at microwave and ~GeV energies ( Finkbeiner|2004 



Dobler & Finkbeincr 2008 Dobler et al. 2010 Su et al. 2010 1 



This scenario requires that the GC has injected an average 
10^® erg/s in CRs over ^ 5 billion years; in this paper we 
set out independent modelling evidence that supports the 
notion that the GC supernovae inject, as required, of order 
this power into non-thermal particles. 



In the second paper of this series ( [Crocker et al.|2010b| 



Paper II) we showed that the GC loses most of the high- 
energy particles it accelerates - precisely as required by 
Paper I. Specifically, we determined the ratios, -Rradio = 
^"yLh/i'SlToh and 7?Tcv = L?.^=v/J^TcvN between the GC's ob- 
served non-thermal radio continuum (RC) and ~TeV 7-ray 
emission and that expected on the basis of the region's FIR 
output (which traces its current star-formation rate) and 
the supposition that the region is calorimetric or a 'thick- 
target' to CRs. These ratios are given by i?radio — 10~^ and 
i?Tev — 10~^ (both uncertain by a factor ~2). On the basis 
of this we found that the GC launches a powerful wind. This 
wind serves to advect most of the CR ions and electrons ac- 
celerated in the GC from the region, launching them into the 
Galaxy-at-large. The identification of this outfiow, however, 
was contingent, upon GC supernovae acting with typical ef- 
ficiency as accelerators of CRs (losing about 10% of the 10^^ 
erg mechanical energy per SN into CRs); in this paper we 
show that, indeed, all the evidence is that GC supernovae 
do act with at least such typical efficiency. 

In this paper we present single-zone, steady-state mod- 
elling of the inner Galaxy. We explain below why both 
the single-zone and steady-state assumptions are physically 
plausible though we emphasise that, with respect to the 
GC's TeV 7-ray emission, others - notably the HESS dis- 
covery paper ( Aharonian et al. 2006 1 - have adopted a 



non-steady-state interpretation. Our modelling allows us to 
render a number of the qualitative conclusions reached in 
our previous papers more quantitative and show how our 
modelling of the GC's particle astrophysics allows us to tie 
down the ranges of a number of parameters including the 
mean gas density and magnetic field experienced by CRs 
and the speed of the outflow that advects them from the 
GC. Our modelling allows us to delimit the potential con- 
tribution made by secondary electrons (and positrons) to 
the region's non-thermal radio emission and determine the 
efficiency with which GC SNRs accelerate CR ions and elec- 
trons. We also discuss at some length the consequence of the 
picture of the GC ISM and its non-thermal particle popu- 
lation that we have developed in the course of the above 
papers and the modelling described below. 



1.1 Structure of paper 

The remainder of this paper is structured as follows: in !|2] 
we delineate the exact region around the GC with which we 
are concerned and describe in detail the RC and ~TeV 7- 
ray emission from the GC we seek to understand. In !|3]we 
set out our detailed modelling of the region's high-energy, 
steady-state particle population and its non-thermal, radia- 
tive output. This affords us an independent handle on the 
region's ISM parameters and confirms the existence of an 
outflow as previously suspected. In ij4] we describe the re- 
sults of our modelling and in i|5] we consider the many im- 
plications of our modelling. These include the determination 
of GC environmental parameters (average gas density, mag- 
netic fleld where CRs are interacting, etc); the result that - 
independent of any uncertainty surrounding supernova rate 
or SNR CR acceleration efficiency - the GC loses ~ 10^^ 
erg/s in non-thermal particles into the Galaxy-at-large; and 
the role this release of CRs plays in explaining large-angular- 
scale non-thermal emission around the GC. In !|6]we present 
our conclusions. Finally, in extensive appendices we describe 
current knowledge about environmental conditions in the 
GC; set out a number of independent determinations of the 
star-formation and consequent supernova rate in the inner 
200 pc; and show that supernovae occur in this region oc- 
cur on average every ~2500 years. We also review the many 
pieces of pre-existing evidence for a powerful GC outflow and 
set out the expectations for its characteristics as informed 
by observations of other star-forming galaxies. 



2 INPUTS 

2.1 Observations of regions of interest 

2.1.1 TeV -ray data: HESS field defined 

The most compact region we concern ourselves with is 
deflned by observations with the H.E.S.S. Imaging Air- 
Cherenkov 7-ray Telescope (hereafter, 'HESS'). After re- 
moval of point TeV sources coincident with Sgr A* and the 
SNR G 0.9-1-0.1 (of angular size ^ 0°.l), this ins trument 
was shown to have detected ( Aharonian et al.|2006 l diffuse, 
~TeV 7-ray emission distributed more-or-less symmetrically 
around the GC over the region defined by \l\ < 0.8° and 
|6| < 0.3° with an intensity of 1.4 x 10"^° cm"^ eV"^ s"^ 



sr^ at 1 TeV (the solid angle denoted by the inner rectangle 
in Fig.[l]Q Only limited and dimmer diffuse TeV emission is 
detected outside this region. Hereafter we denote this region 
as the 'HESS field'. Modelling the particle astrophysical pro- 
cesses that occur within the region seen in projection is the 
chief motivation of this paper. 

Of note is that the spectral index of the GC diffuse 
~TeV emission, at 7 = 2.29 ± O.OT.t^t ± 0.20,^, (where 
the spectrum is defined to be F.y{E-f) oc E~'') is signifi- 
cantly harder than the spectral index of the CR ion pop- 
ulation threading the Galactic disk (i.e., the population we 
detect at Earth) and the HE, diffuse 7-ray emission they gen- 
erate. Disk CRs experience energy-dependent confinement 
(because of the energy-dependence of cosmic-ray scattering 



1 The HESS field-of-vicw is - 
the instrument is 11.3'. 



5° and the resolution (FWHM) of 
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on the Galaxy's turbulent magnetic field structure) and their 
resulting steady-state distribution is, therefore, steepened 
from the injection distribution into the softer ~ ^-2.75 gpgg_ 
trum observed at earth (see, e.g., Aharonian et al.|2006 l. In 
fact, the GC TeV 7-ray spectral index (and that inferred for 
the parent CR ions) is close to that inferred for the injection 
spectrum of Galactic disk CRs, itself within the reasonable 
range of ~ 2.1-2.2 expected (Hillas 2005') for Ist-order Fermi 
acceleration at astrophysical shocks (associated, e.g., with 
supernova remnants). 

Also of significance is the fact that the diffuse TeV 7- 
rays and the molecular gas column (as traced by CS line 
data) are correlated over the central ~ 2°. This supports 
the notion that the 7-rays originate in hadronic collisions 



cal extension of the DNS (with an approximate half- height of 
~ 150 pc) corresponds closely to that of the Galactic centre 



between CR ions and ambient hydrogen (Aharonian et al. 
2006). This correlation, however, breaks down at Galactic 
longitude I = 1°.3. The HESS collaboration (Aharonian et 



al. 2006 1 showed that a very plausible explanation of the 



combined facts of the hardness of the diffuse TeV emis- 
sion and the breakdown in the correlation at / — 1°.3 is 
that the emission is a non-steady state phenomenon ener- 
gised by a single explosive injection of CR ions. This event 
should have occurred too recently for detectable diffusion 
steepening of the CR ion population to have occurred or 
for the leading edge of the diffusion sphere to have reached 
I = 1°.3. Moreover, given the geometry, this event would 
have to have occurred close to the actual GC and, there- 
fore, plausibly in association with the super-massive black 
hole Sgr A* . For a diffusion coefficient typical for that deter- 
mined in the Galactic plane, the extent of the TeV emission 
could be used to infer a time since the explosive injection 
event of ~ 10* years (this transport timescale being much 
less that the synchrotron cooling timescale of the ~ 10 TeV 
electrons required to generate the diffuse TeV emission via 
IC emission, rendering further support to the notion that 
the emission is produced by CR ions). 

There are steady-state alternatives to the general sce- 
nario described above, however, which we will explore in 
this paper. Generically, these alternatives require that there 
be an energy-independent escape or loss process acting on 
the cosmic rays whose characteristic timescale constitutes 
the smallest relevant timescale for the system. Given the al- 
most energy-independent pp loss timescale, collisions with 
ambient gas could satisfy this requirement. Alternatively, 
advective transport affected by a GC outflow or wind might 
define the shortest relevant timescale. 



2.1.2 Radio observations of GC diffuse emission: DNS 
defined 



Radio observations at 74 and 330 MHz ( LaRosa et al 
|2005[ ) have revealed a diffuse but discrete, non-thermal radio 
source ('DNS') approximately covering the elliptical region 
of 3° semi-major axis (along the Galactic plane) and 1° semi- 
minor axis, centred on the Galactic dynamical centre (cor- 
responding to physical scales of ~ 420 pc and ~ 140 pc at 
the assumed GC distance of 8 kpc); refer Fig.[T] T his region 
easily encompasses the CMZ ( [Morris fc Serabyn||1996| , the 
agglomeration of dense molecular matter along the Galactic 
plane within ~ 200 pc of the actual GC and corresponds, 
more or less, to the usual definition of the Galaxy's Nuclear 



lobe (see Sofue & Handa (19841; Law (20101 and references 
therein) and these structures may be more-or-less identical 
and, as is investigated below, due to an outflow from the 
central region of the Galaxy. There is also a suggestive cor- 
respondence between the shape, size, and spectrum of the 
DNS and the extended radio continuum halos seen around 
many star-forming galaxies, e.g. the region of extended, non- 
thermal radio emission surrounding the star-burst galaxies 
NGC 253 ([Mohan et al.|2005||Zirakashvih fc Volk|2006t and 



M82 ( [Seaquist fc Odegard[[199T l 



From an examination of archival radio data we have re- 
cently shown ( [Crocker et al. 2010a I that the spectrum of 
the DNS is non-thermal up to ~10 GHz with an overall 
radio spectral index, oc v~°' of a ~ 0.7 (excluding the 
74 MHz datum which is problematic because interferomet- 
ric and free-free-absorption-affected). A pure power law is a 
poor description of the spectrum (excluded with 3.4a confl- 
dence), however, which in fact exhibits a spectral steepening 
at ~1 GHz (see Fig. [2]| of ~ 0.6. Our analysis of this radio 
spectrum - in concert with ~GeV 7-ray observations by the 
EGRET instrument ( [Hunter et al.[[1997[ ) - has shown that 
the large-scale magnetic field of the DNS region is at least 
50 /xG (radio data is described in detail in the Appendix, 
the analysis of this data is described at length in the Sup- 



plementary Material to Crocker et al. 2010a I. Again using 



archival radio data, we have measured the radio spectrum 
of the HESS field defined above. This is, as for the DNS, 
non-thermal up to ~10 GHz (see Fig.[2| with an overall ra- 
dio spectral index a ~ 0.5 (excluding the 74 MHz datum), 
indicating a very flat synchrotron emitting electron popu- 
lation (and, sub-dominantly, pollution from Hll regions as 
discussed at length below). The spectrum is reasonably de- 
scribed by a power law (again excluding the 74 MHz datum). 

The total radio flux densities measured within the DNS 
and HESS fields receive - in addition to the truly diffuse 
emission we are concerned with - contributions from dis- 
crete sources within each field and are also polluted by 
synchrotron emission from CR electrons along the line-of- 
sight but out of the actual GC. How we deal with these 
back/foregrounds is described at length in section [3] (also 
see Supplementary Material to Crocker et al.[[2010a[ ). 



2.1.3 GeV'j-ray data 

EGRET observations ( Hunter et al.[[T997 l show a pedestal 
of intensity across the inner tens of degrees of the Galactic 
plane. No obvious, diffuse source on the same size scales 
as the DNS radio structure in evident in these data (e.g.. 



Crocker et al. 2010a I and the data, therefore, define only 



Bulge ( Launhardt et al. 2002 Ferriere et al. 2007 \ . The vertl 



upper limits to the GeV intensity that may be predicted by 
modelling of the region. 

In contrast, it is apparent (see Fig. |3| that the Femi 
instrument observes diffuse emission from a region very sim- 
ilar in extent to the DNS. This 7-ray emission is, however, 
again polluted by the contribution of i) individual point 
sources within the field ( [Cohen- Tanugi et al.[[2009| ) includ- 
ing a GeV source positionally coincident with Sgr A and 
the GC TeV point source ( Chernyakova et al.|2010" l and ii) 
by emission from CRs in the line-of-sight along the Galactic 
plane but out of the GC ( Vitale et al.|2009 |. Apparently be- 
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Figure 1. Total intensity image of ttie region at 2.4 GHz. Radio 
map ( [Duncan et al.|1995| convolved to a resolution of 1°. 2x1°. 2 
with contours at 400, 450, 500, 600, 800, 1000, and 1500 Jy/beam. 
There is a striking constancy in the appearance of the radio struc- 
ture from 74 MHz to at least 10 GHz. The large ellipse traces the 
diffuse, non-thermal radio emission region first identified at 74 
and 330 MHz | |LaRosa et al.|["2005[ |. At the assumed distance to 
the GC of 8 kpc, 1° corresponds to a distance of ~140 pc. The 
small rectangle delineates the HESS field. 



Figure 3. Fermi raw count map of emission above 1 GeV around 
the GC (as intensity) with contours showing 5 6 8 10 12 14 18 
counts (14 months' data-taking) smoothed to the Fermi PSF of 
0.3° at 1 GeV illustrating {only) the 7-ray morphology. The ellipse 
and box are again the DNS and HESS fields. The image was made 
using the standard Fermi data reduction pipeline. For the total 
spectrum of the region we employ the data presented by |Meurer| 
|]( |2009^ ; the GeV point source coincident with {l,b) = (0,0) was 
recently described by [Chernyakova et al.| l [2010[ |. 
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Figure 2. Radio spectra for regions under consideration (inten- 
sity vs. frequency). The triangular (blue) data points refer to the 
DNS region (see [Crocker et al.|2010a| and references therein and 
Appendix [dJi . The circular (red) data points refer to the HESS 
field. The vertical dashed line divides the (interferometic) 74 MHz 
data — which do not receive a contribution from the Galactic plane 
synchrotron foreground/ background — from all the (single dish) 
other data (which do). The 74 MHz data are also strongly affected 
by free-free absorption. Note that the 330 MHz flux density for 
the HESS region ( [Law|2010| is indicative only as it suffers from 
an uncertain absolute calibration. It is therefore not used in our 
fitting described below. 



cause of the difficulties attendant upon the removal of these 
back/foregrounds, the Fermi collaboration has not yet pro- 
duced an estimate of the GC's diffuse, GeV 7-ray emission. 
Total GeV intensity spectra have appeared ( Meurer [2009 1, 
however, and these define upper-limits to what our modelling 
may predict. The spectrum of this total emission over the 
inner GC (1° x 1° box centred on the {l,h) = (0,0)) is 



well described by a sharply broken power law that breaks 
downwards from 7 ~ 1.4 to 7 ~ 2.6 at Ecut — 1-6 GeV 



( jMeurer 20091. Such a spectrum is consistent with milli- 
second pulsars (msps) in the GC field (and along the sight) 
dominating the observed emission. The spectrum of indi- 
vidual msps is well-fit as a broken power law with aver- 
age parameters ( |Abdo et al.|[2009) [Malyshev et al 



20101 



7 = 1.5 ± 0.4, Eaut = (2.8 ± 1.9) GeV and L = 10' 
erg/s. The number of such pulsars in the GC environment 
is rather badly constrained given observational difficulties; 
broadly, the Fermi GeV data require ~ few hundred msps 
in the GC field. 



2.1.4 Spatial correlations between GC emission in various 
wavebands 

Although to be the subject of a forthcoming, in-depth study 
by some of u^ we mention in passing here that there ex- 
ists a good spatial correlation between TeV and ~GHz RC 
emission over the GC (also see |Erlykin fc Wolfendale|2007[ ). 
Furthermore, emission in both these wavebands is also cor- 
related with the CS molecular-line data of ITsuboi et al.l 



( [1999 1, particularly after optical depth effects are taken into 
account. These facts support the notion that the observed 
emission is truly from the GC and that the single-zone as- 
sumption we employ (i.e., that CR electrons and protons are 
generating the observed synchrotron and 7-ray emission in 
the same environment) is a justified simplification. 



See Jones, Aharonian, Crocker et al., forthcoming. 
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2.2 Environmental conditions 

Environmental conditions in the GC are reviewed in ap- 
pendix |A] In general the HESS field encompasses a region 
of very high volumetric average molecular gas density (~ 120 
cm~^ as inferred from the data presented by 



Ferriere et al. 



20071, but with most H2 actually concentrated into the low- 



filling-factor cores of the region's giant molecular clouds. 
There is now convincing evidence for a second, relatively low 
density (~ 100 cm~^) but warmer molecular phase envelop- 
ing the denser core material. The dominan t filling-factor 
phase is a very hot plasma of ~ 10'^~^' K (Ferriere et al 



2007 Higdon et al. 2009 and references therein). The re- 



gion's magnetic field is also very strong (> 50 /xG: Crocker 
|et al.||2010a[ ). Finally, we note that this paper and Papers I 
and II reach the conclusion that there is a fast ( ^ 100 km/s) 
outfiow from the GC. There is considerable prior evidence 
for such an outfiow which is summarised in appendix ^C] 



2.3 Timescales: system in (quasi-) steady state 

An important consideration in our modelling is whether the 
system is in steady-state so that, in particular, the typical 
time between injection events (i.e., supernovae) - given the 
supernova rate - can be taken to be negligibly small in com- 
parison with other timescales of relevance, in particular, the 
transport and loss timescales. In fact, given the supernova 
rate determined in Appendix[B] this condition is met as sum- 
marized in Fig. |4] where the various displayed timescales can 
be (approximateljj^ summarized as 
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(Here the ionization and bremsstrahlung loss time have been 
adjusted to account for heavier nuclei with Z > 1 in the 
background gas.) Note the hadronic (pp) loss time is only 
slightly longer than the bremsstrahlung loss time with the 
same direct dependence on uh- Ionization losses become im- 
portant for low energy protons but, in the energy range of 
interest to us, pp is by far the dominant loss process. The 
above establishes that the system meets the requirement for 
being in (at least) quasi-steady state (so that we may seek 
steady-state solutions of Eq[2] below). 
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Figure 4. HESS region timescales for i) (horizontal solid band) 
the inverse of the supernova rate; ii) (dashed horizontal lines) par- 
ticle escape (assumed to be advective with energy-independent 
velocity: lines are, from top to bottom, for cases of v^i„^ = 100, 
300, 1000 km/s); iii) (red lines) electron cooling (for cases of B = 
10, 100, 1000 ^lG) in a) = 10 cm-^ (solid); and b) volu- 
metric average nn — 120 cm"^ (dotted); iv) (blue lines) proton 
cooling in a) njj = 10 cm"'' (solid); and b) volumetric average 
riH — 120 cm~3 (dotted). Note that proton energy-loss timescales 
are always longer than electron energy-loss timescales over the 
pertinent particle energy range and that calorimetry generically 
requires tjoss < tesc- Particle energy loss processes are described 



2.3.1 A note on the TeV "/-ray deficit at I — 1°.3 

As described above, the HESS collaboration interpreted the 
break-down in the correlation between the diffuse TeV emis- 
sion and the molecular gas column density as indicating a 
non-steady-state situation, viz., the 7-ray emission had been 
energised by a single CR ion injection event close to the ac- 
tual GC, the leading edge of which had not yet reached this 
longitude. Within the context of the steady-state modelling 
we perform this explanation is not available to us: this deficit 
must speak of either a spatial change in CR ion population 
and/or gas characteristics at this distance from the GC. In 
fact, such is not unreasonable: molecular studies (Tanaka et 
|al.||2007] ) have revealed a distinct 'I = 1°.3 complex' which 
is located (in projection) at the geometrically-privileged po- 
sition defined by the eastern terminus of the CMZ. This 
complex contains gas characterised by unusually large ve- 
locity extents and a 'puffy' structure ( Oka et al.|1996[ ) with 
at least nine expanding shells, some possessing enormous 
(~ 10(51-52) gj.g-) kinetic energy and displaying associated 
SiO clumps. The exact interpretation of this phenomenol- 
ogy is a matter of debate: the structure might have been 
energised by a very recent burst of massive star-formation 
(I Tanaka et al. 20071. Alternatively, the I = 1°.3 gas may 



For full expressions see Crocker et al. I 2007a I 



have only recently arrived in the central region and collided 
with the central condensation ( [Sawada et al.|[2004[ ). In any 
case, given these anomalies, a break-down in the correlation 
between TeV emission and gas density at ? = 1°.3 - while 
certainly interesting - does not necessarily reveal that a sin- 
gle event energised the region's CR ion population; steady- 
state interpretations of the 7-ray phenomenology and the 
particle astrophysics behind it are possible (the simplest 
probably being that the I — 1°.3 gas complex gas is not 
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even contiguous with the rest of the CMZ, but only appears 
so because of projection effects; we will examine these issues 
at greater length in future work). 



3 MODELLING OF NON-THERMAL 
PARTICLE POPULATIONS 

We have created a one-zone model of the injection, cooling, 
and escape of relativistic protons, electrons, and secondary 
electrons (and positrons) from the HESS region. As justified 
above, in our modelling we assume that the particle astro- 
physics can be accurately described to be in quasi-steady 
state and that the particle transport timescale is energy- 
independent. 



3.1 Modelling 

3.1.1 Steady-state particle distributions 

We assume homogeneity and isotropy in our one zone mod- 
elling. In this case the transport equation is; 



dnJyEx 
dt~ 



d 



Ex nx{Ex 

Tia.s{Ex) 



nx{Ex] 
r.s.(Ex 



+ QJyEx 



(2) 



(Ginzburg & Syrovatskii 19641, where Qe[Ee) denotes the 
injection rate of particles of type x, x £ {e,p} and we ac- 
count for both particle escape and energy loss over, respec- 
tively, timescales of Tesc and Tioag. The steady-state distribu- 
tion of non-thermal particles can be written then as 



rio.s(Sx) 
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E' 



dE'x TiassjEx 

E'J r^E'^) 



For our purposes a sufficiently accurate approximation to 
the above can be written: 



nx{Ex 



TiasB{Ex)T^, 



- Qx{Ex 



(3) 



where 7 is the spectral index of the (assumed) power-law (in 
momentum) proton or electron spectrum at injection, and 
we assume the escape time, t^sc, to be energy-independent 
(see discussion in Appendix [e||. Here we take it that the ac- 
celeration timescale is less than both the loss and escape 
timescales up to some large _E"'"' (otherwise the steady- 
state cosmic ray proton distribution would not likely be 
a power law, in contradiction with the expectation sup- 
plied by the 7-ray data). Note that Eq. [s] reduces to the 
usual leaky-box expression, nx{Ex) — r^BcQxiEx), in the 
limit that t^^^ <^ tiobb and to the 'thick target' expression, 
TLxjEx) ^ Tia^,/{-^ — l) Qx{Ex), in the opposite hmit (see, e.g., 
[Crocker et al.|[2007a[ ). Given the energy-independent adia- 
batic energy loss time and the almost energy-independent pp 
cross-section, Eq. [S] renders a very accurate approximation 
to the steady-state proton distribution. For electrons, over 
the parameter space of interest to us (especially for spectral 
indices of concern which are 2.0 or steeper) we find from 
numerical evaluation that Eq. [3] reproduces Eq. [3] to better 
than ~ 35%. 



3.1.2 Cooling Processes 

Relativistic protons and electrons may both suffer adiabatic 
energy losses in doing PdV work on the escaping wind fluid 
(with the requisite momentum exchange presumably com- 
municated dominantly via the CRs' scattering on Alfvenic 
disturbances in the magnetized wind plasma: e.g., |Socrates| 
et al.|2008 l which are given by: 

[Ecr) = Ecr~~ \t,sc EcR . (4) 

rfi 3 3 

In addition, as already alluded to, there are cooling processes 
that affect relativistic electrons or protons exclusively (see 
[2.3 1. For protons, note that the steady-state distribution - 



given the near energy-independence of the pp cross-section 
over the apposite energy range (and the assumed energy- 
independence of the escape time) -will have a spectral shape 
identical to the injection distribution. The inelastic collisions 
suffered by the protons lead to the production of neutral 
pions (which decay into 7-rays) and charged pious whose 
decay products include relativistic electrons and positrons. 
Self-consistently, our model also tracks the radiative output 
of these "secondary electrons" . 

In the case of electrons (and positron^ we have that 



rL.iE.)^-E./^(E.) 



(5) 



where the total loss rate is given by the sum of all loss pro- 
cesses acting to cool the relativistic electron distribution, 
viz. adiabatic deceleration, ionization, bremsstrahlung, and 
synchrotron and inverse-Compton emission: 



dE, 
dt 



{E, 



EdK 
dt 



(6) 



where i £ {inztn, brems, synch, IC}. Each of these processes 
has a characteristic and different energy dependence as de- 
scribed in §2.3| Roughly, the cooling rate is independent 
of energy in the case of ionization, oc E, in the case of 
bremsstrahlung, and oc E^ in the case of synchrotron and 
inverse Compton emission: full expressions can be found in 
[Crocker et al.| ( |2007a[ ) . These different energy dependences 
introduce spectral features into the steady-state electron dis- 
tribution that are in addition to the features inherited from 
the injection distribution and occur at an energy where one 
cooling process takes over from another (or from escape). 
Note, however, that over the range of energy over which ei- 
ther bremsstrahlung or (energy-independent) escape is the 
dominant process, the spectral index of the steady-state elec- 
tron distribution is the same as that at injection. 



3.1.3 Radiative processes 

The same cooling processes described above lead to the ra- 
diation of photons that are potentially observable (though 
it does not follow, of course, that, e.g., for electrons syn- 
chrotron radiating at ~GHz frequencies these particles are 
only or even dominantly cooled by synchrotron emission). 



Note that the ultra-rclativistic energy scale of the leptons under 
consideration implies that energy losses for positrons and elec- 
trons can be taken to be identical with high accuracy. Positron 



annihilation can also be ignored: see, e.g., Crocker et al. (2007a 
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We self-consistently calculate the radiation of the relativis- 
tic (primary and secondary) electron population via syn- 
chrotron, bremsstrahlung, and inverse Compton emission. 
We also calculate the 7-ray emission induced by the hadronic 
collisions experienced by the relativistic protons. Expression 
for these processes are given by Crocker et al. (2007a I. Note 
that, in calculating the synchrotron cooling and emission, 
we assume an isotropic electron distribution. This implies 
the solid-angle-averaged magnetic field perpendicular to the 
electron's direction of travel is {B±} = n/4 ~ 0.78 B. 



of the emitting electron) synchrotron is not a relevant pro- 
duction process for ~TeV 7-rays in this environment. Also 
note that we will ignore 77 attenuation which in principle is 
non-negligible at energies of ^ 100 TeV for the photon-rich 
GC environment ( [Zhang et al.|2006|poskalenko et al.|2006 \ 
but which introduces an attenuation which is small on the 
scale of other uncertainties in our modelling. For pp 7-ray 
production, F^^y{E-,), we use an interpolation between (at 
< 100 GeV) the results of 



E. 

E-y > 100 GeV) the resuhs of 



Kamae et al. 



Kelner et al. 



20061 



( 2006 1 and (at 



3.1.4 Signal at radio wavelengths 

The (intrinsic, unabsorbed) signal we wish to model at radio 
wavelengths consists solely of the synchrotron emission due 
to primary (le) and secondary electrons (2e): 

Fsynch{l^) = Flynchi^) + F^ynchi^) (7) 

where the flux densities of these separate contributions are 
functions of the following parameters (defined previously) 
over which we minimize the function set out below: 

Flynch{v) = F}y„^,,{ly,Qe,"/s,tesc,B,nH) 



Faynchi^) — F^ynchiv, Q2e.,tesc, B , U h) 



(8) 



where Qe(TeV) 



„Qp(TeV), Q2e = Q2e{Qp,1p,nH), and 



we assume the injection spectral indices obey 7e = 7p = 7. 
Note that, for clarity, here and below we do not list other 
in-principle parameters that are fixed in our modelling (dis- 
tance to, volume of, and solid angle of source region). Also 
note that at scaling energies 720 — 7p but we do not assume 
this relation. Rather, the spectrum of secondary electrons 
(resulting from the collisions of a particular, power-law pop- 
ulation of CR protons and taking into account the kinemat- 
ics of the TT — >■ /i — )■ e decay chain) is as given by a 
MATHEMATICA interpolation of particle yields from ac- 
celerator data on charge d pion production in pp collisions 
(as described at length in Crocker et al.|[2007a[ p| 



3.1.5 Signal at ^-ray energies 

At TeV energies we have to consider emission by both 
(highly) relativistic electrons and protons. Relevant pro- 
cesses are neutral meson (mostly tv'^) decay following 
pp collisions for protons and inverse-Compton (IC) and 
bremsstrahlung for both primary and secondary electrons. 



FreviE^) = F^^y{E^) + F^ly{E^) + F}ly{E^ 



(9) 



where 



F^.v{E^) = Fi" 

+ FtcV 



,nH) 



esc 1 I 

Tev(-^7' Qe, 7^ ^esc, -S, Uh) 
■brems 



{E^,Qe,'y,tesc,B,nH) 



(10) 



Again, note that the interstellar radiation field - an im- 
portant parameter in describing IC cooling and emission - 
is fixed in our modelling with an energy density of ~ 90 
eV cm""^ and spectrum as described above. Over reason- 
able parameter values (in particular, the requisite energy 



Data kindly provided by Todor Stanev. 



3.2 Back/Foregrounds at radio wavelengths 

As well as the diffuse, synchrotron emission that we seek to 
determine from our modelling of the region's total RC out- 
put, the measured fiux density also receives a significant con- 
tribution from i) the foreground/background of synchrotron 
radiation by relativistic electrons in the Galactic plane along 
the line of sight (but out of the GC), the so-called Galac- 
tic Synchrotron Background (GSB; note that the 74 MHz 
datum is supplied by the The Very Large Arrajj^an inter- 
ferometer, which is not sensitive to the large-angular-scale 
GSB at this frequency); ii) discrete radio sources within the 
field; and iii) free-free (thermal bremsstrahlung) emission 
from Hll regions in the region and in the Galactic plane in 
front of the GC. 



3.2.1 Galactic plane synchrotron background 

The Galactic plane contribution to the total radio emission 
is significant - possibly approaching half the measured emis- 
sion - so we must deal with this component. We normalize 
our GSB model at 408 MHz via recourse to the Galactic 



plane survey of Haslam et al. ( 1982 1 and, following the pro- 



cedure of LaRosa et al. (20051, choose a constant longitude 



slice free of discrete sources (well outside of the DNS) near 
I = 354°. 5 to measure the latitude dependence of the 408 
MHz flux density normalization. At higher frequencies we 
expect the GSB to have a spectral index of a = — 0.70±0.12 
( Platania et al.||2003[ ) with I„ ex u". Note that the GSB is 
expected to display a break at a frequency of ~ 22 GHz 
(V61k l989l, well above the region where we are fltting. As 
we cannot directly measure the GSB at higher frequencies 
due to the presence of other backgrounds, we leave the GSB 
spectral index as a free parameter in our fitting. This compo- 
nent of the detected emission is described, then, by a single 
fltting parameter: 

FGSB{y,lGSB) = Fgsb(408 MHz)(^ ^^g ^^J .(11) 



3.2.2 Free- free emission and absorption 

Given the compact and extended Hll regions within the GC 
and along the line-of-sight ( LaRosa et al. [2005 1 we expect a 
contribution from free-free (thermal bremsstrahlung) emis- 
sion to the total radio flux. Self-consistently, we must take 

^ The VLA is operated by the National Radio Astronomy Obser- 
vatory (NRAO). The NRAO os a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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into account the optical depth due to free-free absorption 
in modelling the detected radio flux due to synchrotron and 
free-free emission at the source. Our modelling of free-free 
processes works in the following way: at sufficiently high fre- 
quency, the free-free optical depth, r// becomes small and 
free-free emission, F^^ , might account for anything up to 
100% of a high-frequency flux datum. Our procedure, 
then, normalizes the total free-free emission at 10.29 GHz 
where Tfj<^l for reasonable parameters. At this high fre- 
quency, the free-free optical depth and emission - connected 
through the emission measure - are related via: 



r// (10.29 GHz) 



F^-'' (10.29 GHz) 
STRESS 5(10.29 GHz,r) 



(12) 



where we make the simplifying assumption that the Hll re- 
gion covers the entire HESS region solid angleQ f^HESS, and 
-8(10.29 GHz, T) is the intensity of a black-body of tem- 
perature T at 10.29 GHz. We adopt T = 5000 K in this 
calculation as motivated by radio recombination line obser- 
vations of the GC ( [Ferriere et "aL][2007| [Law et al.||2009| ). 
Self-consistently, the free-free optical depth at any frequency 
can then be determined as 



rffH = r„(10.29 GHz 



)( — 
' V 10.29 



GHz 



(13) 



and, assuming the geometry of an intervening screen, the 
detected free-free emission at any frequency is then 



Fffiv) = Ohess 5(10.29 GHz,r) (l-e-""*"') 



(14) 



On the other hand, the synchrotron emission region - lo- 
cated behind the assumed free-free screen - has its flux at- 
tenuated according to 



(15) 



where F, 





ynch 



is the flux that would be detected from the 



region in the absence of free-free absorption. 



3.2.3 Discrete sources in the field 

We do not attempt to subtract off the radio flux associated 
with discrete sources in the fleld (whi ch we have previously 



Crocker et ar]|2010a 



determined to be at the ^ 10% level 

consistent with the expectation that individual SNR con- 
tribute 10% of the total radio flux from normal galaxies: 
Condon 1992, Zirakashvili & Volk 2006) because the TeV 
7-data over the given region does not resolve out any such 
sources except for the Sgr A point source (and the radio 
flux associated with the Sgr A source - in contrast to its 
TeV flux - is small in comparison to the region's integrated 
radio flux; e.g., 20 Jansky at 1.4 GHz, cf. the almost 2000 
Jy over the entire region: Crocker et al.|2007b l. 



Note that, in fact. Ho absorption over the field is somewhat 
patchy; [LaRosa et al.| l [2005[ l. Thougii not strictly justified given 
we find (dof < 1 already for best-fit regions of the current pa- 
rameter space, we have investigated the introduction of a further 
parameter for the Hll covering fraction in our modelling. We find 
that a covering fraction 60% is required for a good fit; this is 
consistent with the results of jMezger &: Pauls| l |1979[ l who found 
that most thermal RC emission from the region was contributed 
by extended, low-density ionized gas (discrete Hll regions con- 
tributing the remainder). 



3.3 Modelling procedure 

Our approach is to find - as a function of environmental and 
other parameters - the steady-state populations of relativis- 
tic protons and electrons within the region of interest. (Note 
that in our modelling we neglect for simplicity the poorly- 
constrained ionic component of the CR hadronic population 
heavier than protons.) We then self-consistently determine 
(for the same environmental parameters) the radiative out- 
put of these populations, as described above. Finally, we use 
a x'^ minimization procedure to determine the parameters 
describing the proton and electron populations whose radia- 
tive outputs in the given ISM environment give the best fit 
to the total radio flux (density) and the ~TeV 7-ray flux 
detected from the HESS fleld (as previously noted, the GeV 
data only deflne upper limits to the broadband emission 
from the region and we do not attempt to fit them directly) . 

Environmental parameters that vary within our mod- 
elling are magnetic field _B, ambient hydrogen number den- 
sity riH (in whatever form), and the energy- independent 
timescale over which particles are advected from the sys- 
tem, t^sc = h/vwind, where h = 8 kpc tan(0.3°) ~ 40 pc. 

We assume that protons and electrons are injected into 
the GC ISM with distributions governed by power laws in 
momentum (with identical spectral indices 7p = 7e = 7)- 
The relative normalization of the injection distribution of 
electrons to that of protons (at relativistic energies) is given 
by the coefficient Kep which is also left as a free parameter 
in our modelling. The absolute normalization of the distri- 
bution of protons at injection is, finally, specified by Qp (in 



units eV cm 



, also a free parameter within the model 



3.3.1 Chi- squared function 

We determine the volume of parameter space for which 
an adequate fit to the data is obtained via a numerical 
minimization within MATHEMATICA (using the internal 
FindMinimum function, with method Interior- Point) of a 
function: 



X = xidio + Xtcv + Xgsb 



where 

2 ^ V- (F„dio(!^i)°'' - -Frad.o(l^O)' 

Xradio 



(16) 



(17) 



with 



{u,} = {74 MHz, 1.408 GHz, 2.417 GHz, 

2.695 GHz, 8.35 GHz, 10.29 GHz} 

and 

F..d.o(i^>) = e-^"'^''\Fl^„,n{y) + F^^^^^iu)) 

+ Fgsb{u) + F;j[u) (18) 

and where 



XtoV 



with 



(19) 



FTev(i?^) = F^„v(5^) +F^UE,) + F^UE^) . (20) 
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The parameters over which we minimize the x ^re 
Qp, K,ep,'-f,v^ind, B,nH, and 'jgsb- Guided by considera- 
tions set out in !|C] where necessary to slef-consistently pa- 
rameterize the wind, we assume a thermalization efficiency 
in the range 0.1 ^ rj ^ 0.9 in our modelling. 



4 RESULTS OF MODELLING OF HESS 
REGION 

4.1 Typical parameters 

Our fitting procedure finds acceptable fits to the data we 
model for magnetic field amplitudes around ~ 100 ^G, wind 
speeds around «„ind ~ fewxlOO km/s, total power in all 

^39 



10 erg/s, gas densities around 



^, injection spectral indices 7 ~ 2.4, ionization 
and electron to proton ratios (at 



non-thermal particles of 
uh ~ 10 cm 
rates of ^ 

injection at 1 TeV) of k^p ~ lO"'^. (The implications and the 
statistically acceptable range of each of these are discussed 
at greater length below.) 

Below we illustrate (figs. [5] and [6| the modelled broad- 
band and radio spectra for the best-fit case (which achieves 
~ 7.9 for dof = 9). We also show (Fig.^ the inferred CR 
proton and elctron spectra - at injection and steady-state - 
for this particular case. 



4.2 Pure 'hadronic', pure 'leptonic' scenarios 
excluded 

As part of our modelling we have examined whether the 
broad-band, diffuse emission from the HESS region might 
be explained either within pure 'hadronic' or 'leptonic' sce- 
narios. In the former case, synchrotron emission by sec- 
ondary electrons and positrons, by construction, explains 
most of the non-thermal RC emission while the non-thermal 
TeV emission is mostly accounted for by neutral pion decay 
(with a subdominant contribution from secondary electron 
IC and bremsstrahlung) . Note that Thompson et al. (20061 
have claimed an important - even dominant - role for sec- 
ondary electrons in generating the RC detected from star- 
bursting systems. In a 'leptonic scenario', primary electron 
synchrotron explains the RC while IC and bremsstrahlung 
by higher-energy members of the same primary population 
explains the TeV signal. In fact we find poor fits for both 
these limiting cases; the best-fit pure hadronic case has 
X^/dof = 14.9/7, acceptable only at ~ 2.1a and requires 
a very strong (~4 mG) field; the best-fit pure leptonic case 
has X^ /dof — 28.3/7, acceptable only at ~ 3.7a and requires 
a very weak (~16 /iG) field and a very tenuous (~ 2 cm~'^) 
gas environment. 



4.3 Existence of nuclear wind confirmed by 
modelling 

Pure modelling of the non-thermal spectrum does not con- 
strain the wind speed but the imposition of physically- 
motivated constraints (as described in the Fig. [s] caption) 
on the parameter space shows that the wind speed is > 200 
km/s and most likely several hundred km/s (see Fig. |8|. 
Given the 1.4 x lO**" erg/s mechanical power supplied by 




Figure 5. HESS field broadband spectrum energy distribution 
for the best-fit case. Fitted parameters include magnetic field am- 
plitude (^ 200 ;iG), ambient gas density (6 cm~^), at-injection 
eloctron-to-proton ratio at 1 TeV (~0.005), and wind veloc- 
ity (~400 km/s). Curves are divided into (i) dashed - pri- 
mary electron emission; (ii) dotted - secondary electron (and 
positron) emission; and (iii) solid — total emission at a given 
photon energy. Emission processes are: blue - synclirotron; red 
- bremsstrahlung; green - inverse Compton; and brown, dot- 
dashed - neutral meson decay. The total 7-ray flux is shown in 
purpl e. Data are from (at ~ IQ-^) radio, (at B-^ ~ 2 keV) 
Ginga jYamauchi et al.|l990| with , (at ~ 1 GeV = lO'' eV) 
Fermi | |Meurer |2009|l, and (at Ej ~ 1 TeV = lO^^ eV) HESS ob- 
servations l | Aharonian et al. |2006l l . We only display modelled syn- 
chrotron emission from CR electrons at radio wavelengths; other 
processes that combine with this to give the region's total radio 
emission are shown in Fig. [6] (so the overall flt to the radio data is 
better than it appears here). Also note that the Ginga- measured 
X-ray flux is, apparently, dominantly due to optically-thin ther- 
mal plasma emission of temperature ~ 10^ K. The Ginga Large 
Area Proportion Counter FOV is well-matched to observations 
on the size scale of the TeV emission and, in fact, show that the 
putative thermal plasma covers a very similar solid angle to the 
HESS region. Finally, note that the Fermi data points are for to- 
tal emission from the region and deflne only upper limits to the 
diffuse emission we seek to model. 



supernovae, a conservative (because assuming 100% ther- 
malization and no mass loading) upper limit on the wind 
speed is ~ 1200 km/s. 

4.4 GC magnetic field determined by modelling 



As remarked upon at length in Appendix |A2[ the GC mag- 
netic field has been uncertain by about two orders of mag- 
nitude for some time though our recent work ( [Crocker et al.| 
[2010a| | has determined a lower limit to the field on the size 
scale of the DNS of 50 /iG. Relaxing all constraints and fit- 
ting only to the broadband spectrum we do not significantly 



constrain this range. If we impose the reasonable (e.g., Volk 
et al.|1990 1 - but admittedly heuristic - constraint that the 



ratio of proton to electron differential number density in- 
jected into the ISM satisfy K^p < 0.01, we find that the 
magnetic field is required to satisfy B > 50 fiG at 2a 
confidence with a minimum at ~160 /iG. Requiring that 
the outflow velocity is governed by Eq. |C1| in addition to 
other, reasonable constraints (see caption of Fig|9|, tends 
to pick out best-fit B values around this value. It is inter- 
esting to remark, moreover, that imposition of these extra 
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Figure 6. Radio data and the modelled radio emission due to 
various process for the best-fit parameters described in the text. 
Note that the 74 MHz datum does not receive a contribution from 
line-of-sight Galactic plane synchrotron (GSB) emission, hence 
the (unphysical) step in the modelled total emission curve (solid 
blue curve). Other curves are:- solid, purple: total synchrotron; 
dotted, purple: total synchrotron in absence of free-free absorp- 
tion; dashed, black: primary electron synchrotron; dotted, red: 
secondary electron synchrotron; dotted, green: GSB; and dot- 
dashed, brown: free-free emission. 



EE 



(a) Modelled proton and elec- 
tron populations at injection 
(both assumed to be governed 
by power laws in momentum) 
for the best-fit case. 



(b) In situ charged Icpton 
spectra (as formed by the 
loss and escape processes de- 
scribed in the text) for the 
best-fit case. The rising spec- 
trum of secondary electrons at 
very low (and radiatively unde- 
tectable) kinetic energies is due 
to 'knock-on' electrons. The lo- 
cal electron spectrum is as de- 
scribed by [Aharonian et al.| 
l |2009| |. 



(c) In situ proton spectrum (as 
formed by the loss and escape 
processes described in the text) 
for the best-fit case. The lo- 
cal proton spectrum is as de- 
scribed by |Webber| | |1998| |. 

Figure 7. Inferred particle spectra pertaining to the best-fit case. 




2.0 

Log[v„i„d/km/s] 

Figure 8. -)(^ jdo^ as a function of v^i„d- Solid curve is for a 
obtained from a pure fit to the broadband data as described 
in the text {dof = 9). The dashed (purple) curve follows the 
introduction of a term into the foi' power injected into CRs 
(central datum value set to power expected on basis of SN rate 
plus stellar winds, 1.4 X 10^^ erg/s) and constrains the power 
advected from the system by CRs to be enough to generate the 
DNS emission, E^dv > 10^* erg/s {dof = 10). The dotted (yel- 
low) curve follows the imposition of the heuristic (e.g., |Volk et al.] 
|1990[ l constraint that the ratio of proton to electron differential 
number density injected into the ISM satisfy Kep < 0.01 and also 
that njj ^ = 6 cm~^, where n^-i- is the volumetric aver- 

age over the HESS region of all gas in plasma phase as inferred 
from the data of'Ferriere et al.| ( |2007[ l. The dot-dash (green) curve 
shows a fit to the broadband data which also attempts to repro- 
duce -Rtov — 0.01 and -Rradio — 0.1 (see Paper II and jjl]) with 
riH > ~ 6 cm"''. The vertical dashed and dot-dashed lines 

show, respectively the approximate maximum (~ 1200 km/s) and 
minimum wind speeds according to Eq. |C1| with thermalization 
efficiencies of »7j^^J^ = 1 and i?™"^,,, = 0.1. The vertical dotted 
line shows the approximate gravitational escape speed from the 
mner GC of ~ 900 km/s | |Muno et aL,2 004) and the vertical solid 
line shows the mid-plane sound speed, Ca = ^ 7 Pc/ pc, in the 
wind plasma for central values of E and M. 



constraints requires the field be ^ 300 /xG, less than the 
~600 /.iG necessary for the magnetic field alone to ensure 



hydrostatic equilibrium (see Appendix A2.2 1 



This seems to tell against a |Thompson et ar] ( |2006| l-type 
scenario, at least as it might apply to the GC: the magnetic 
field in the GC - while certainly much larger than that in 
the disk - is neither large enough to ensure the region is 
calorimetric to electrons nor (alone) to establish hydrostatic 
equilibrium (though one should keep in mind that the region 
is not in hydrostatic equilibrium in any case given the evi- 
dence for an outfiow which is driven by the combined pres- 
sures of the plasma and non-thermal components). If the 
situation really is as described by [Thompson et aL] ( |2006[ ) in 
dense star-bursts, this is, then, a point of difference between 
the GC and these systems. 

The field is also significantly weaker than the ~ mG 



1987 



amplitude suggested by analysis ( Yusef-Zadeh et al 

[Morris fc Yusef-Zadeh(T989t of the region's non-thermal fil 
aments (NTFs; see Appendix A2|. 
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-4.0 -3.5 
Log[B/Gaus,s] 

Figure 9. ^ function of B. The solid curve is for a 

obtained from a pure fit to the broadband data as described in 
the text (dof = 9) with wind speed as an independent variable. 
The short dashed curve is a pure fit to the broadband data as 
described in the text {dof = 9) with wind speed as an independent 
variable and Kep < 0.01. The dotted curve requires the wind 
speed obey Eq. |Cl| with »7therm ^ fr^e parameter (constrained to lie 
in the range 0.1 < »?therm < 0.9) with the constraint that the total 
non-thermal pressure (generated by the modelled magnetic field 
and CR population) be no larger than the pressure launching the 
wind as given by Eq. |C3| (cio/ = 9). The dot-dashed curve also 
requires the wind speed obey Eq. |Cl| with r^tjie^m ^ free parameter 
and, in addition, introduces a term into the for power injected 
into CRs (central datum value set to power expected on basis 
of SN rate plus stellar winds, 1.4 X 10^^ erg/s) and constrains 
the power advected from the system by CRs to be enough to 
generate the DNS emission, Eadv > 10'^* erg/s (dof = 10). The 
long dashed curve constrains the magnetic field and CRs to be 
in energy-density equiparition under the assumption that both 
these phases are volume-filling. The vertical dotted lines show the 
magnetic field amplitude i) in equipartition with the interstellar 
radiation field (TSRF'), ii) in equipartition with the very hot 
plasma ('plasma'), iii) required for electron calorimetry (defined 



by r 
cm" 



and conservatively assuming njj = { nu)yai ~ 120 
760 km/s as supplied by Eq. CI and adopting 



"cool — '-esc 
-3 yfriin 

wind 

'^tho™ = 0-f)i ^'^'f iv) required to achieve byitself hydrostatic 



equilibrium. 



4.5 Power demanded by modelling 

From direct modelling of the broadband spectrum we find no 
clear minimum in the x^y expressed as a function of power 
injected into all CRs; see Fig. [lO] This figure, however, also 
shows the effect of introducing the dual constraints of i) re- 
quiring that the pressure contributed by CRs and magnetic 



fields (Eq. C3 1 is no larger than the total pressure launch- 



ing the wind and ii) self-consistently setting the wind speed 
to that given by the mechanical power and mass injected 
into the system (while still allowing the thermalization ef- 
ficiency to vary between 0.1 and 0.9 - see Eq. CI I. With 



the imposition of these constraints the allowed parameter 
space is significantly reduced, favouring a power into non- 
thermal particles ^ 10^® erg/s, the value predicted by the 
supernova rate and adopting a fiducial 10^" erg per SN into 
CRs. In other words, we reach the interesting conclusion 
that GC SNRs are at least as efficient as disk SNRs as CR 
accelerators, losing ^ 10% of their mechanical power into 
acceleration of such particles. 
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Figure 10. x^ ^ function of power injected into all CRs and 
CR electrons. The solid curve shows a fit to the broadband data 
only (dof = 9) for total CR power. The dotted curve shows a 
fit to the broadband data only (dof = 9) for total CR electron 
power. The solid vertical lines show the power in non-thermal 
particles assuming they capture a factor rjQji of 0.01, 0.1, and 1 
(as marked) of the total mechanical power released into the HESS 
region by supernovae (and stellar winds), adopting the central 
value of supernova rate of 0.04/century. (The vertical dashed lines 
indicate, for tiqh = 0.1 the range of the power given by the factor 
of ^ 2 uncertainty in the supernova rate.) The dashed curve shows 
a fit to the total CR power constrained such that the wind speed 
obey Eq. |C1| with fjthcrm a free parameter (constrained to lie in 
the range 0.1 < ??thorm < 0.9) with the constraint that the total 
non-thermal pressure (generated by the modelled magnetic field 
and CR population) be no larger than the pressure launching the 
wind as given by Eq. |C3| (dof = 9). The dot-dash curve shows 
the equivalent constraints for the case of CR electron power. 



4.6 Medium where emission is occurring 

On the basis of our modelling we can infer the character- 
istic gas density for the region in which the modelled non- 
thermal radiation is generated. Fitting only to the broad- 
band data we find a very shallow minimum for 1 cm~^ ^ 
nH > 100"^ cm"^ with an even larger nn range acceptable 



at 2a level (see Fig. 111. Note that the volumetric average 
uh in the region is ~ 120 cm~^ with most of the mass con- 
tributed by molecular gas found at very high number density 
( ^ 10* cm~^) in the cores of the region's giant molecular 
clouds. It is clear that the hypothesis that emission origi- 
nates dominantly in these regions is disfavoured. We also 
indicate in the figure the volumetric-average number den- 
sity of gas in all plasma phases, 71^+ , (including the 'very 
hot', 'hot', and, dominantly, 'warm' H''' phases identified by 
Ferriere et al. 2007 see Appendix A1.2|. Given that, ac- 
cording to the analysis of Higdon et al. (20091, relativistic 



particles will not be trapped in the very hot phase (for which 
Uh 55 0.1 cm~^) njj+ arguably constitutes a lower limit to 
the effective, path-averaged nn the region presents to CRs. 

Inclusion of data in addition to the broadband spec- 
trum constrains the allowed uh parameter space in a very 
interesting fashion. The dashed and dotted curves in Fig. 
|11| shows the effect of including terms for iiTov(— 0.01) 
and 7?radio(— 0.1) in the x^ (see Paper II and ^ and re- 
placing the previously-independent wind speed parameter 
with Eq. |C1| with r;t™„ = 0.4 and r/uic™ = 0.9, respec- 
tively {dof = 11) with the maximum thermalization see- 
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Figure 11. ~x? /dof as a function of njj (strictly for dof = 9). 
Solid curve is for a obtained from a pure fit to the broadband 
data as described in tlie text {dof = 9) with wind speed as an 
independent variable. The dashed curve fits to the broadband 
data and also attempts to reproduce Rtov — 0.01 and -Rradio — 
0.1 (see Paper II and with wind speed given by Eg. |C1| with 
0.1 ^ ^therm ^ 1-0 {dof — 10). The vertical dashed line shows 
the nfj+ , the volumetric average rifj over the HESS region due 
to gas in all plasma phases. Th e vertical dot-dashe d line shows 
the density, n^"^* n i r.m-3 



0.1 cm~'=' I Muno et al 
very hot phase plasma phase The vertical do' 



20041, of only the 
ted line shows the 



volumetric average njf over the region for all matter as inferred 
from [Ferriere et al.| | |2007| l . 



nario favouring higher densities and the minimum lower. 
Over ?7J55e?m < '7thorm < Vthc^m the 2o--acceptable region now 
extends from 0.5 cm^'^ < uh < 80 cm^"^, with even the 
volumetric uh value formally excluded at the 2a level. If 
we accept as a plausible lower limit, we find an nn 

in the range 6-80 cm~^ is favoured by our analysis. This - 
together with the best-fit magnetic field described above - 
is consistent with the region of parameter space favoured by 
the heuristic arguments presented in Paper II. 

This result hints that CRs - even those at super- TeV en- 
ergies - do not penetrate freely into dense molecular cloud 



material; cf. theoretical (Dogiel & Sharov 1990 
aL]|2007[) and empirical studies ( [Protheroe et al." 



Gabici et 



20081 IFu- 



jita et al. 2009) of this issua Given the observational fact 



of the angular correlation between both RC and TeV 7- 
ray emission and the molecular material (as traced by CS), 
these results are consistent with the observed radiation be- 
ing generated in low density Hi or Hll material surrounding 
the clouds (cf. |Higdon et al.[ 2009) or, indeed, their relatively 
tenuous and warm molecular envelopes of the cores (but not 
dominantly in the hottest, lowest-density, volume-filling Hll 
phase: cf. Fig. [TT|) . 

These results are also potentially consistent with recent 



findings from Fermi observations of the LMG ( Knodlseder et 
al.|2009[ ). Here it was demonstrated that neutral and molec- 
ular hydrogen templates poorly fit the GeV data whereas an 
ionized hydrogen template (tracing current star formation) 
provided a relatively good fit. In the GC, our analysis sug- 
gests that the synchrotron emission from ~GeV electrons 



and 7-ray emission from super- TeV protons may mirror this 
result. In part, this can be simply explained as a result of the 
fast advection timescale: the non-thermal particles are ad- 
vected from the system in a timescale ~ 4 x 10* t;wind/(1000 
km/s) yr to be contrasted with the advection timescale into 
the giant molecular cloud complexes of tadv = Rcioud/o-v 
which for, e.g., Sgr B {Rdoud = 12 pc and cr„ = 40 km/s) 
evaluates to fadv = 3 x 10^ yr (of course, this comparison 
neglects the poorly constrained diffusion timescale into the 
clouds; [Protheroe et al.|2008 1 

We note in passing that, as stated above, our modelling 
makes the simplifying assumption of a single zone. This im- 
plies, in particular, that the synchrotron emission we model 
- due to ~ GeV energy electrons - can be treated as though 
it is emanating from the same environment as the TeV signal 
which is (predominantly) generated by impacts of multi-TeV 
protons on ambient gas. There is the possibility, of course, 
that this assumption does not accurately refiect the circum- 
stances in the GC; that, in fact, multi-TeV protons sample, 
on average, denser gas than the OR electrons. We have inves- 
tigated this possibility in modelling which allows for different 
environmental parameters {B and uh) in describing the CR 
electron and proton environments. This modelling actually 
prefers, if anything, electron environments slightly denser 
than the proton ones. We take this as a good indication 
that the single-zone modelling described above (fortuitously 
or not) renders a good description of the GC's high-energy 
astrophysics. 



4.7 Injection ratio of electrons to protons 

Fig |12| shows that fits to the broadband spectrum of the 
region only poorly constrain the at-injection ratio of elec- 
trons to protons (at 1 TeV) so that, for instance, good fits 
are possible for either of the limiting cases where, on the 
one hand, equal numbers of CR electrons and protons are 
injected and, on the other hand, equal power is injected in 
CR electrons and protons. As for other parameters, how- 
ever, the imposition of additional, physically-motivated con- 
straints (as described in the figure caption) restrict the pa- 
rameter space somewhat. It is certainly of note that the 



minimized for k^v in the range predicted ( |Bell|1978| ) for the 
case that the injection spectra for both electrons and pro- 
tons are pure power-laws-in-momentum (and equal numbers 
of protons and electrons at thermal energies are injected 
into the acceleration process) with spectral indices in the 
range 7 = 2.2 — 2.3 (consistent with the observed TeV spec- 
tral index); admittedly, however, the 2cr-allowed parameter 
space of this parameter is still fairly large. Also consonant 



10 



with the usual picture (e.g., Thompson et al.|2006 l. Fig 
shows that CR electrons take around 10% of the total CR 
power for the best-fit region of parameter space 



DISCUSSION: IMPLICATIONS OF 
MODELLING RESULTS 



* Also see Jones, Crocker et al. 2010, 'Australia Telescope Com- 
pact Array Radio Continuum 1384 and 2368 MHz Observations 
of Sgr B', submitted to the Astronomical Journal. 



5.1 Cosmic ray heating and ionization of GC gas 

On the basis of our modelling, we find that the CRs in the 
HESS region induce a local ionization rate of ^ 10~^^ s~^. 
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et al. 



indeed seem to be in a two-temperature distribution (Oka 
[T9981 [Rodriguez-Fernandez et al.][200T| [Goto et al 



2008 



Oka et al. 2005 



Yusef-Zadeh et al. 20071 with warm 



low-density gas enveloping cooler, denser cores. 

The exclusion of CRs from the densest molecular core 
regions suggests CRs do not significantly heat or ionize the 
ambient, dense gas where stars actually form. This is an in- 
teresting conclusion to reach in the context that the GC 
be considered a true star-burst analogue; it has recently 
been argued that CR heating and ionization of molecular 
gas in star-burst environments radically alters the condi- 
tions for star-formation (away from those encountered in 
more quiescent star-forming regions) thereby leading to a 



Los 



Figure 12. ^ jdoj as a function of Kep, the at-injection ratio of 
electrons to protons at _E = 1 TeV (dof = 9). The solid curve 
is for a obtained from a pure fit to the broadband data as 
described in the text (dof = 9). The dashed curve requires the 
wind speed obey Eq. |Cl| with r^therm a free parameter (constrained 
to lie in the range 0.1 < »?therm < 0.9) with the constraint that the 
total non-thermal pressure (generated by the modelled magnetic 
field and CR population) be no larger than the pressure launching 
the wind as given by Eq. |C3|(do/ = 9 ). The vertical dashed Unes 
shows the ratio predicted ( |Bell|1978[ l when the injection spectra 
for both electrons and protons are pure power-laws-in-momentum 
with spectral indices 7 = 2.2 and 2.3 (consistent with the observed 
TeV spectral index). 



This is precisely in the range required by the data discussed 
in Appendix |A1| 

Moreover, on the basis of the results obtained by 
[Suchkov et al.| ( |1993[ ) for star-bursting environments (see 
their Fig. 1), this ionization rate would keep the GC molec- 
ular gas at a temperature of ^ 50 K, again precisely in 
the range required by the data (also see [Yusef-Zadeh et ah] 



top-heavy IMF (cf. Papadopoulos 2010 Papadopoulos et 



al. 2010 1. This mechanism would seem to be excluded in 



the GC. Whether it might still operate in 'true' star-bursts 
(which, e.g., are dense enough that calorimetry is putatively 
achieved: Thompson et al.][2006 l we cannot yet say. We do 
note, however, that turbulent velocities typical for the most 



luminous starbursts (which may exceed 200 km/s: Downes & 
Solomon|1998 l may imply advection timescales into molecu- 



lar clouds similar to t^ind in these systems (cf. the GC case, 
S 4.6 1 , potentially implying CRs in these clouds do reach into 
the dense gas (cf. Contini fc Goldman|2010 1. 



5.2 GeV emission from the GC 

It is clear from the broad-band spectrum presented in Fig. [5| 
that we explain ^ 10% of the observed 7-ray emission at 
~GeV energies. This is to be expected given, as remarked 
above, that we only have access to total flux data from Fermi 
and that these receive a substantial contribution from dis- 
crete sources in the GC ( [Chernyakova et al. 12010) and from 
both discrete and diffuse emission sources along the line of 
sight. 



2007 Dogiel et al. 2009b I. On the other hand, given this 



5.3 Explaining DNS/GC lobe radio emission 



temperature, the kinetic pressure in the GC region would 
then imply that the molecular material would typically be 
at ~ 10^ cm"'', a number density higher than that found (or The observational study by Law ( 2010 1 revealed that the 



allowed) by our fitting (cf. fig 111, which is typically much 
more tenuous. 

A compelling resolution of these facts is that the molec- 
ular gas exists, in fact, in two phases that are in pressure 
equilibrium at ~ 5 x 10^ K cm~^: dense cores of 10^ cm""' 
and relatively tenuous envelope material. Note that in this 
scenario pressure equilibrium is established between the ki- 
netic pressure of the core material and the virial pressure im- 
plied by the turbulent motions of the tenuous gas. Also note 
that for a typical GC gas turbulent velocity dispersion of 20- 
40 km/s, pressure equilibrium implies the tenuous material 
has a density in the range 50 — 200 cm"'', comfortably within 
the range given by our modelling. Finally, both phases will 
be heated by CR ionization to similar kinetic temperatures 
assuming the low-energy heating/ionizing cosmic rays can 
actually penetrate into the cloud cores. As we have already 
seen, however, such penetration is extremely unlikely. This 
fact, together with the higher cooling efficiency of the dense 
core material (see fig 1 of Suchkov et al.|1993 1 , suggests that 
it should have a lower kinetic temperature than the envelope 
material. This is all precisely in line with the phenomenol- 
ogy discussed in Appendix |A1[ the GC molecular gas does 



~GHz spectral index of the GC lobe steepens with increas- 
ing distance (both north and south) of the Galactic plane. 
This constitutes strong evidence for synchrotron ageing of a 
CR electron population transported out of the GC, entirely 
consistent with the wind scenario we have been postulat- 
ing here. Also of note is that the vertical extension of the 
GC lobe is very similar to the DNS. We previously postu- 
lated (Paper II), therefore, that the DNS and GC lobe are 
essentially identical structures synchrotron-illuminated by 
CR electrons advected from the inner GC (essentially the 
HESS region) on a wind (cf. [Seaquist fc Odegard[[199T on 
the radio halo of M82). 

Three basic conditions must be satisfied in order that 
electrons escaping from the HESS field go on to pervade the 
DNS region and generate most of its synchrotron emission 
at ~GHz wavelengths, viz.: i) the spectrum of the electrons 
leaving the HESS region (i.e., the steady-state distribution 
in that region as given by Eq[3| must match the spectrum 
at injection required for the DNS electrons; ii) the power in 
electrons leaving the HESS region must be enough to sup- 
port the DNS electron population; and iii) the time to trans- 
port electrons over the extent of the DNS must be less than 
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or of the order of the loss time over the same scale. As shown 
in Paper II, all these conditions are well-met within our sce- 



jects inhabiting or processes occurring within the HESS re- 
gion do, indeed, accelerate ^ 10^^ erg/s in cosmic rays (cf. 



nario. Firstly, previous modelling by some of us ( Crocker et 
al. 2010a I shows that electrons are injected into the DNS 



with a spectral index in the range 2.0-2.4, perfectly consis- 
tent with the spectrum of electrons leaving the HESS region. 
Secondly, the same modelling demonstrated the total power 
required to sustain the steady-state electron population of 
the DNS is ~ fewxlO"^^ erg/s, again consistent with the 
power advected away by HESS region electrons. Thirdly, we 
have also shown previously that the magnetic field of the 
DNS region has a (2 o) lower limit at 50 ^.G and a prob- 
able value of ~ 100 /iG. Requiring that electrons can be 
transported over the entire vertical extent of the DNS (~ 1° 
from the plane, corresponding to 140 pc) implies a wind 
speed of a least a few hundred km/s (see Fig. 1 of Paper 
II), again consistent with the wind speed predicted by the 
modelling of the HESS region. More detailed modelling of 
electron transport out of the central HESS region and the at- 
tendant synchrotron cooling in order to match the detailed, 
position-dependent radio spectrum of the lobe (cf. the stud- 
ies by Zirakashvili fc Volk| |2006) and [Heesen et al.| ( [2009 1 
of NGC 253) we leave to future work. 



5.4 Explaining the WMAP 'haze' and Fermi 
'bubbles' 

From detailed morphological and spectral analysis of large 
angular scale maps of microwave and ~GeV 7-ray emission 
produced, respectively, by the WMAP and Fermi instru- 
ments, Finkbeiner an d co-workers (|Finkbeiner|2004[ |Dobler| 
& Finkbeinerl 120081 iDobler et all 120101 iSu et al.l 120101) 



have made the interesting claim that there is evidence for 
anomalous, non-thermal 'hazes' of emission centred on the 
GC. This emission is apparently also coincident with ex- 



tended soft X-ray emission traced by ROSAT (Snowden et. 
|al.|1997[ ). The most recent analysis ( |Su et al.|2010[ ) shows the 
anomalous ~GeV emission to resemble two bilateral bubbles 
of emission centered on the core of the Galaxy and extending 
to around 10 kpc above and below the Galactic plane. 



As has been explored at length elsewhere (Crocker & 
Aharonian]|2010| ), the WMAP haze and the Fermi bubbles 
are naturally explained within our model as due to the pop- 
ulation of relic cosmic rays protons and heavier ions (not 
CR electrons as has mostly been explored in the litera- 
ture), injected in association with the high areal density star- 
formation in the Galactic center. In this scenario, the sec- 
ondary electrons created by the hadronic collisions of these 
ions - occurring within the rarefied wind environment - nat- 
urally supply the synchrotron emission putatively detected 
by WMAP. The accompanying neutral meson decay (with a 
sub-dominant contribution from the IC emission by the sec- 
ondary electrons), on the other hand, naturally explains the 
Fermi emission. This explanation requires that these cosmic 
ray hadrons i) carry off a power of at least 10^^ erg/s and 
ii) have been injected over an extremely long timescale, at 
least 5 Gyr (the latter in order that pp energy losses in the 
extremely rarefied outfiow, uh ~ 0.01 cm~^ put the local 
CR ion population into saturation). 

An important result of this current work - in support 
of the explanation of the large-angular-scale microwave and 
GeV emission we advocate - is the demonstration that ob- 



Fig. 10 1. We have made the case here that GC supernovae 
can supply the power required to sustain the acceleration of 
this non-thermal particle population but note that this de- 
termination (that ^ 10^^ erg/s is injected into CRs) is in- 
dependent of uncertainties surrounding the supernova rate, 
efficiency of SNRs as cosmic ray accelerators in the GC en- 
vironment, etc. Note, in fact, that we cannot exclude, in 
principle, the possibility that other sources in the region 
might ultimately supply the power required by the CRs; 
the super-massive black hole Sgr A* is an obvious candi- 
date ( Aharonian et al~]|2006 1 particularly as recent analysis 
( Chernyakova et al.|2010 1 of the GeV and TeV point sources 
coincident with it suggest that it, too, roughly supplies of 
order 10^^ erg/s in cosmic rays. 

The long timescale of injection required to explain the 
haze/bubbles would - were it associated with any other 
Galactic source of structure - immediately call this explana- 
tion into question. The morphology of the observed bubbles 
(Su et al. 2010 1, however, unambiguously associates them 



with the GC, probably the only structure in the Galaxy for 
which spatially-localized star formation over such a long pe- 
riod is credible ( Serabyn fc Morris||1996 Figer et al.|2004 1 . 

Another interesting aspect of this scenario is that the 
wind of plasma and CRs leaving the GC should neither es- 
cape to infinity nor fall back to the Galactic plane. In fact, 
given a wind speed of few hundred km/s, it is probable that 
the wind injected flow out of the GC stalls at a height of a 
few kpc - just as required - as has recently been investigated 



at length by Rodriguez-Gonzalez et al. ( 2009 see Appendix 
Cil. 



5.5 Galactic bulge positrons 

A nuclear wind could also help disperse positrons into the 
Galactic bulge (cf.|Dermer fc SkibojiggTllCasse et al.|2004| 
Prantz os 2006; ITotani||2006| [Higdon et al.||2009| IGebauerl 
& de Boer 2009: Churazov et al.|20T0 l, potentially explain- 



ing the morphology of the observed 511 keV emission line 



( Knodlseder et al.|2005 Weidenspointner et al.|2008 l as we 
also explain in Paper I. It should be noted, however, that 
this postulate does not address the problem of where the 
positrons come from in the first place. 

A prime candidate process to supply the requisite, low- 
energy positrons is inverse beta decay of radioisotopes in- 
jected by GC Type la supernovae. It has been claimed, how- 
ever, that GC SNla can, only supply 10% of the positrons 
required by observations ( [Schanne et al.|[2007[ ). This deter- 
mination assumes 3.3% of the positrons from ^^Ni decay 
escape the Type la envelope. This fraction, however, ob- 
tained from normalizing to a statistical sample of optical 
light curves ( Milne et al.|1999 l, is model-dependent and de- 
pends, in particular, on assumptions about the magnetic 
field strength and topology and the state of the gas in the 



expanding remnant (Milne et al. 1999 Higdon et al. 2009 



Martin et al. 2010). Indeed, higher values for it have been 



countenanced (Higdon ct al.__2009) and direct upper limits 
are at the ^ 12% level for a number of disk SNRs (Martin 



et al. 2010). In any case, the highly porous and, in parts, 



extremely rarefied GC environment may allow many more 
positrons to escape the envelopes of GC SN la's than typical 
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for disk-type environments (cf. Parizot et al.|20d5 Higdon et 
al.|2009 Martin et al.|20l"0 i. This consideration may extend 
to other classes of SN in the GC environment. 

The many massive, young stars in the GC also consti- 
tute a compelling source population for an unusual sub-class 
of SN Type Ic, thought to originate in Wolf-Rayet explo- 
sions and exemplified by GRB-associated SN 2003dh, that 
have been speculated (Cass e et al.|2004 l to be sources of the 
bulge positrons, as have GRBs more generally (Bertone et 
al.|2006 Parizot et al.|2005 1 . Another possibility is that pro- 
cesses unique to the GC may be acting to supply positrons 



et al. 20061. Also of note is the fact that particle trans- 



(e.g., tidal disruptions of massive stars by the SMBH: Cheng 
et al.||2006| |Totani"200 6|) as, inde ed, may the population of 
millisecond pulsars ( ^Wang et al.||2006 l that the GeV spec- 
trum of the region indicates may be present. We leave in- 
vestigation of these speculations to further work. 



As a final note in this section, Churazov et al. (20101 



have recently concluded (on the basis of a spectral study 
of SPI/INTEGRAL data on the circum-GC 511 keV line 
emission) that the annihilating positrons are injected into 
a hot (~ 10^ K), volume-filling ISM plasma phase which 
undergoes radiative cooling. The initial temperature of this 
gas is somewhat cooler than favoured by our analysis but 
we do note that, broadly, this scenario would seem to be 
consonant with our picture. It is also of interest that the 



prior modelling by Rodriguez-Gonzalez et al. ( 2009 1 of a 



star-formation-driven GC outflow with initial temperature 
~ 10^ K shows that the plasma cools to ~ 10^ K by the 
time it has reached a height of ~ kpc. According to the 



modelling of Churazov et al. (20101 at such a temperature 



the positrons may annihilate so, again, the combination of 
these two ideas would seem to point to a model that matches 
the observed ~kpc extension of the 511 emission. 



6 CONCLUSIONS 

We have modelled in detail the high-energy astrophysics of 
the inner ~200 pc of the Galaxy with the aim of explaining 
a number of observations including the region's non-thermal 
radio and ~TeV 7-ray emission. On the basis of our mod- 
elling/analysis we determine the following:- 

(i) Supernovae occur in the inner 1°.6 (in diameter) of 
the Galaxy at a rate of approximately 1 per 2500 years and 
release ~ 1.3 x 10*° erg/s in mechanical power (which is sup- 
plemented by massive stellar winds at the ~10% level; see 



Appendix Bl I. This is precisely enough to drive the fast out- 
flow ( ^ 100 km/s) from the region and support the region's 
non-thermal particle population. Given the data, moreover, 
are consistent with the current star-formation and super- 
nova rates being close to their long-term average values, a 
system in quasi-steady-state is suggested. 

(ii) Despite the GC resembling a star-burst in many ways. 



a number of the circumstances Thompson et al. ( 2006 1 claim 
for star-bursts apparently do not pertain in the GC: 

(a) We conform that, in general, 'calorimetry' fails in 
the GC; il, 95 % of the power injected into non-thermal 
particles is advected from the system, viz. ~ 10^'' erg/s. 
This despite the high surface mass density of gas and the 



port is overwhelmingly advective; given, presumably, the 
strong, turbulent magnetic fields of the region, particle 
diffusion timescale are much longer than either advection 
or energy-loss timescales for both protons and electrons. 
This severely limits the role of the super-massive black 
hole on larger angular scales consistent with the deter- 
mination that supernovae and stellar winds can drive the 
high-energy activity (setting aside the 7-ray point source 
coincident with Sgr A*). 

(b) Despite the high volumetric-average gas density, 
secondary electrons probably play only a minor role in 
generating the observed radio emission from the GC. Both 
pure leptonic and pure hadronic scenarios for explaining 
the broadband emission of the HESS region are rejected 
with reasonable confidence, mirroring the conclusions of 



Crocker et al. ( 2007a I on the Sagittarius B region, one of 



the largest giant molecular cloud complexes in the CMZ. 

(c) The large-scale magnetic field of the GC is bounded 
by 60 ^iC ^ B > 400 /iG. While certainly strong with 
respect to the typical Disk field amplitude and, in fact, 
comparable to the field strengths determined for the cen- 
tral regions of nearby starbursts (e.g., NGC 253: Heesen 



et al. 20091, this is short of the strength required, by it- 



self, to establish either electron calorimetry or hydrostatic 
equilibrium (though, as has been noted, hydrostatic equi- 
librium does not seem to pertain in any case given the 
existence of the wind - driven by the combined pressure 
due to all ISM components). 

In this context, one should keep in mind that the extent to 
which Thompson et al.'s starburst model pertains in reality 
to genuine starbursts is contested. The recent TeV 7-ray de- 
tection of NGC253, e.g., has been interpreted by the HESS 
collaboration, on the one hand, as indicating this system 
falls well short of being calorimetric ([The HESS Collabora-| 
tion: F. Acero|2009 l but Lacki et al. (20101 reach a different 
conclusion (though the latter authors also grant that the 
system does fall short of being fully calorimetric). We also 
note with interest a very recent paper by |Lacki fc Thomp-| 
(20101 investigating the high-energy astrophysics of a 



number of starburst-like environments including the GC (in 
the context of evaluating the contribution of electron syn- 
chrotron to the diffuse, soft X-ray emission detected from 
these systems; the authors find a negligible contribution in 
the GC conclusion with which we agree: cf. Fig.[5f. 

(iii) The GC magnetic field is also significantly weaker 
than the ~mG amplitude predicted by analysis of the NTFs 
( Yusef-Zadeh et al.|[T987| [Morris fc Yusef-Zadeh||1989| ). On 
the other hand, it is likely very close to energy equipar- 
tition/pressure equilibrium with the very hot plasma (and 
turbulent motions of the dense molecular gas) . The best-fit 
magnetic field also falls very close to the value required to es- 



tablish the phenomenologically-expected (Thompson et al 



2006 1 scaling Ub — 2 Uisrf- These facts tend to suggest 



;.g., a 



relatively high magnetic field amplitude (cf. Thompson 



the field is generated locally in the GC ISM through, e.] 
turbulent dynamo. 

(iv) The modelling tends to favour regions of parameter 
space where the wind is driven approximately equally by the 
pressures ascribable to thermal (plasma) and non-thermal 
(magnetic field and CRs) GC ISM phases (cf. |Zirakashvili fc| 
Volk|2006[ on NGC 253). The wind has an effective density 
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of nu ~ 0.03 cm~ , a characteristic speed near the mid 
plane satisfying 100-f 200 km/s and most probably 400-800 
km/s and advects about 0.02-0.03 Mo/year of gas from the 
GC (see Appendix B5l. 

(v) Given the total power implied by our modelling and 
the region's SN rate, the simplest - though somewhat sur- 
prising - global explanation of the GC's non-thermal emis- 
sion is that the region's SNRs are the dominant source of cos- 
mic rays which they accelerate with at least typical efficiency. 
This despite the much larger volumetric-average higher- 
density ambient medium which would be naively expected to 



significantly alter the evolution of the SNRs (Spergel & Blitz 



1992 ). This is consistent with the picture that the GC super- 



novae are exploding into a highly porous medium ( Reynolds 
et al. 2008 Ksenofontov et al.|2010 1 excavated in part by the 



progenitor star's wind ( Volk fc Biermann|1988 Biermann & 
Cassinelli|1993|l, stellar winds in general, a nd/or by previous 
supenovae (cf. [Westmoquette et al.|2009 l. 

On the other hand, it should be borne in mind that the ev- 
idence cissociating GC SNRs with acceleration of cosmic rays 
is only circumstantial and, in fact, any process proportional 
to the supernova rate and producing an average ~ 10^*^ erg 
per SN in non-thermal particles satisfies the phenomenol- 
ogy (cf. Volk et al.||1989 19901. Thus, it remains possible 
that at least some acceleration takes place, e.g., on general 
GC ISM turbulence injected by the collective activities of 
supernovae and stellar winds ( Drury et al.|200T l. 

It is also interesting that - given the play afforded by ob- 
servational uncertainty - we have no strong evidence from 
the overall energetics and the FIR luminosity for (or against) 
any strong biasing of the IMF of the GC stellar popula- 



tion (cf. Persic & Rephaeli 2010 Papadopoulos 2010 Pa- 
padopoulos et al.|2010 |. 

(vi) CR electrons accelerated in the HESS region lose 

10 % of their power to synchrotron emission there, generat- 
ing thereby the diffuse, non-thermal radio emission detected 
across the same region. 

(vii) Protons lose 1 % of their power to pp collisions 
on ambient gas in the inner few degrees. These collisions 
explain the diffuse, TeV emission detected across the Central 
Molecular Zone (electrons only making a minor contribution 
to this emission). Interestingly, these pp interactions occur 
in a relatively low density gas phase ( ^ 30 cm~'^) despite 
the fact that most ambient gas finds itself at much high 
density ( ^ 10'* cm~"^) in the cores of giant molecular clouds. 
This is consistent with the picture that diffusive transport 
timescales in the GC (in this case, into the clouds) are much 
longer than advective ones. 

(viii) The Coulomb collisions of the low energy members 
of the in situ, cosmic ray population (both protons and elec- 
trons) can inject the heat to maintain the temperature of 
the anomalously hot, relatively low-density molecular gas 
and supply much of the unusually high ionization power re- 
quired by various observations. On the other hand the molec- 
ular cores are at a lower kinetic temperature probably as a 
combination of their higher cooling efficiency and, presum- 
ably, the fact that even high energy CRs do not penetrate 
them during the time they remain in the region. This seems 
to rule out a role for CR-modification to the region's dense 
gas chemistry (cf. Papadopoulos|[2010 Papadopoulos et al 
2010l ). 

(ix) With plausible physical constraints, our modelling 



favours (at 2g level) CR proton energy densities in the ap- 
proximate range 2.5 eV cm~^ < C/p < 50 eV cm""^ with a 
best-fit value around 20 eV cm~^. The allowed range of this 
parameter reflects the a priori uncertainty in the typical uh 
value where CR ion collisions are taking place. The lower end 
of the acceptable Up parameter range - corresponding to the 
upper end of the allowed gas density (i.e., tih — >■ {nH)vo\) 
- is only a factor of a few larger than the local CR energy 
density, and matches that found in the original analysis of 
the HESS group ( Aharonian et al. 2006 1 where it was as- 



sumed that all of the molecular gas was 'accessible' to the 
TeV CRs. 

Likewise, our modelling favours (at 2a level) CR electron 
energy densities in the approximate range 0.2 eV cm~^ < 
(7e < 4 eV cm"''. The range of this parameter reflects the 
a priori uncertainty in the magnetic field amplitude where 
the synchrotron emission is occurring. 

The existence of the fast wind is what ensures that the 
in situ CR density at the GC is probably only a few times 
larger than the solar value at Ecr — GeV - despite an areal 
density of SNRs (and massive, windy stars) in the region 
much larger than typical for the Galaxy at large - and also 
explains why the (non-thermal) radio continuum emission 
from the GC is smaller than expected on the basis of the 
region's total FIR luminosity and the FRC (Paper II). 



(x) We find that 'equipartition' ( Beck & Krause 2005 ) be- 
tween the magnetic field and the non-thermal particle pop- 
ulation can be satisfied over at least some of the parameter 
space favoured by our modelling. Formally, good fits consis- 
tent with equipartition are possible up to above 10~* Gauss 
or magnetic field and CR proton energy densities surpassing 
250 eV cm~^. Such fits, require, however very small ambient 
gas densities (1 cm~'') and vanishing wind speed: these seem 
quite implausible to us in the light of all other evidence. We 
think it likely that, in fact, equipartition does not hold in 
the GC ISM. 

(xi) The relativistic electrons that escape the region on 
the wind can explain the large-scale, diffuse radio emission 
detected from the Galaxy's nuclear bulge {\l\ < 3°), in par- 
ticular their synchrotron emission explains the non-thermal 
spectrum of what has been labelled in previous work the 
DNS ([LaRosa et al.|2005[ ) and the GC lobe ( |Sofue fc Handa 
p84l|Law|20T0f ; 

(xii) We have not explained the majority of the 
apparently-extended ~ GeV emission from the region but 
much of this emission is due to point sources certainly in- 



cluding a source coincident with Sgr A* ( Chernyakova et al 



20101 with preliminary indications that, not unexpectedly. 



millisecond pulsars also make a large contribution to the 
extended emission. 

(xiii) If - as consistent with the work of others based on 
the current stellar population ( Serabyn fc Morris|1996 Figer 
|et al.|[2"004| [FTgcr 2008) or upper limits to the pulsar pop- 
ulation ( Lazio fc Cordes||2008 1 - the current SFR has been 
sustained for multi-Gyr timescales, the wind emerging from 
the GC and the CR ion population it carries compellingly ex- 
plain the large-angular-scale, non-thermal emission detected 
at microwave and ~GeV energies. 

In summary, the following picture emerges from the 
studies described here and in the preceding Papers I and 11: 
the gradual accretion of gas into the Galactic nucleus - pre- 
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sumably under the action of the Galactic stellar bar, possibly 
assisted by dynamical friction and magnetic viscosity effects 



dBinney et al.'1991';'Bl itz et al.|19 93"Serabyn fc Morris'1996| 
Morris fc Serabyn,1996||Beck et al.^l999,,Stark et al._2004P 



over timescales approaching the age of the Galaxy - has led 
to sustained star-formation inside a region circumscribed by 
the Galaxy's inner Lindblad resonance. Inside this radius lie 
the vast (but dynamic) gas reservoir represented by the CMZ 
and the nuclear stellar population resulting from the persis- 
tent star-formation. Concomitant with this star-formation 
activity are supernovae, whose remnants accelerate the GC 
cosmic ray population and heat the region's hot plasma; the 
tell-tale, diffuse TeV 7-rays and X-rays produced by these 
two extend across the same inner region of the Galaxy. 

Another by-product of the star-formation is the large- 
scale super- wind that emerges from the GC. This acts to 
keep the energy density of the various GC ISM components 
in check, in particular, advecting CRs to very large scales 
above and below the Galactic plane. Acting to also advect 
the magnetic field, the wind may well also have a crucial 
role in forming the NTFs ( [Shore fc LaRosa|1999|[Boldyrev fc| 
Yusef-Zadeh 2006) ) which, intriguingly, are also only detected 
over the same inner region of the Galaxy as the phenomena 
referred to above. In the picture of ( [Shore fc LaRosa|1999[ ), 
the NTFs do not trace a global poloidal field topology, but 
instead are dynamic structures, akin to cometary plasma 
tails, formed from the interaction of the large-scale, mag- 
netized plasma outflow draping the region's dense molecu- 
lar clouds. Consonant with the above, the latest NIR po- 
larimetry measurements ( [Nishiyama et al.|2010[ ) apparently 
reveal a global field topology that is toroidal in the plane 
(for \l\ < 2.0°, \b\ < 0.4°, closely corresponding the HESS 
region), bending to almost perpendicular for |6| > 0.4°, pre- 
sumably in the wind zone. 

Finally, although it should be remembered that, of 
course, some stochastic variation in the background star- 
formation rate is expected, we have pointed out here that the 
current star-formation rate in the GC seems to be close to 
its long-term average value. This apparent stability suggests 
the operation of self-regulation or feedback mechanisms in 
which many of the actors discussed here - the magnetic field, 
the CR population, the wind - undoubtedly play some part, 
but we leave the detailed working-out of this to further work. 
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APPENDIX A: THE GC ISM 
Al Gas distribution 

The interstellar medium throughout the GC region consists 
of a number of co-existing phases with very different physical 
parameters which we now partially describe. For full reviews 
- particularly of the Hi and relatively cooler plasma phases 
not described below - the reader is referred to the work of 



Ferriere et al. (20071 and Higdon et al. (20091. 



A 1.1 Molecular gas 

In its molecular phase, the inner ~ 500 pc of the Galaxy 
contains ~5% of the Galaxy's allocation of molecular hy- 



drogen (a total mass of a few times lO'^ Mq: Ferriere et al 



2007), mostly within the asymmetric condensation along the 



Galactic plane known as the Central Molecular Zone (CMZ; 
Morris fc Serabyn||1996[ ). The overall geometry of the CMZ 
structure - seen in projection - is very difficult to recon- 
struct. Much of the dense molecular gas seems to lie on a 
~180 pc radius ring aroxmd the GC with a relatively tenuous 
interior (see Morris fc Serabyn|1996 Launhardt et al.|2002 



Liszt|2009 and references therein), akin to the nuclear star- 



forming rings commonly observed in barred spiral galaxies 
nominally similar to our own. Inside this radius a number of 
giant molecular cloud complexes are to be found, including 



Sgr B2, the most massive in the Galaxy (see Crocker et al, 
2007a and references therein). 



Aside from these morphological issues, the CMZ is char- 
acterised by: 

• Unusually high inter- and intra-molecular clump veloc- 
ity dispersions of 15-30 km/s ( Oka et al.|19 98 Gusten et al. 



|2004[ ), ranging to even higher values > 50 km/s in partic- 



ularly active regions like Sgr B2 (see, e.g., Crocker et al 



2007a and references therein); cf. typical disk cloud veloc- 
ity dispersions of ~ 5 km/s. 

• Unusually high molecular densities: a typical H2 
molec ule will find itself in an ambient density of ~ 10'* 
( Paglione et al.|[l998 l; cf. typical disk cloud densities 



of ~ 10 cm~ - though it should be emphasised that the 
filling factors of the extremel y dense molecula r cloud cores 
are extremely small, 0.01 (Oka et al. 20051, and the gas 



distribution is likely highly fragmented (rather akin to the 
situation within luminous star-bursts: e.g., |Westmoquette et| 
ar][2009t . 

• Unusually high kinetic temperatures (molecular mate- 
rial is typically found with temperatures in the range 50-70 
K but co-existing phases with temperatures of up to 300 K 
are detected; cf. typical disk cloud temperatures of ~ 15 K) 
and virial pressures (of ~ 5 x 10® K cm~^ - cf. that typical 



for the Galactic disk, ~ 2 x 10" K cm ; [Spergel fc Blitz 



1992) 



Pervasive shocks, as revealed by thermal SiO emission 



( Martin-Pintado et al. 1997 Huettemeister et al. 1998 ) and 



as revealed by surveys of molecular lines and dust continuum 



emission at mm and sub- mm wavelengths ( Oka et aL]|1998[ 
Tsuboi fc Miyazaki|[T998l [Pierce-Price et al.||2000| |Oka et| 



al. 2010 1 - filaments, arcs, and shells, indicative of local. 



turbulent sources ( Tanaka et al.|2007) and explosive events. 



In addition to the dense molecular cores, in the last few 
years evidence for another, comparatively diffuse (~ 100 



cm ^) and hot (~ 250 K) molecular phase has emerged 
from com parison of the inte nsity of different e mission lines 
of *^C0 (|Oka et al.||l998| ), other molecules ('Rodn'guezJ 
Fernandez et al.|2001 1, and from spectroscopic observations 



of line absorption by ( Goto et al.|2008 Yusef-Zadeh et 
al.|2007 Oka et al.|2005 1 at infrared wavelengths. This envl 



ronment, which represents ~30% of the total molecular gas 
by mass ( Ferriere et al.p007 l, appears to be unique within 
the Galaxy to the GC and requires a distributed heating 
mechanism. The large columns observed towards the 
GC furthermore imply an ionization rate of Si lO"'^^ s~*, 
much higher than typical in the diffuse ISM of the Galactic 
disk (~ 3 X 10"^'^ s"*). These same observations, moreover, 
show - contrary to the expectation from chemical model cal- 
culations - that the relative number density of is much 
higher in diffuse clouds than in dense clouds towards the 
GC. In contrast to the gas temperature, dust temperatures 
throughout the CMZ are anomalously low at around 20 K 
( Pierce-Price et aL||2000 ), suggesting that a mechanism for 
gas-heating is required for the CMZ unrelated to ultra-violet 
radiation. 



A1.2 Hot plasma and X-ray observations 

In its plasma phase. X-ray continuum and Fe line obser- 
vations apparently reveal a two-temperature plasma con- 
taining 'hot' (~ 1 keV) and (mysteriously) 'very hot' (6- 
9 keV) components ( Koyama et al. [1989 Yamauchi et al. 



fT990l p^aiieda et al.|[l997| |Muno et ai.|20d4[ [Belanger etld 
,2004[ ). The X-ray emission from the putative very hot com- 
ponent is strongly concentrated within the inner ~150 pc 
(in diameter) of the Galaxy ( Yamauchi et al.|1990 Belmont 



et al. 



2005[ ). As first observed by [Spergel fc Blitz| ( |1992^ 
there may be pressure equilibrium (at (3 — 6) x 10" K cm~"^: 
Koyama et al.[[T996 Muno et al.[[2004 l between the kinetic 



pressure of the putative very hot plasma phase and the virial 
pressure implied by the turbulent motions of the molecular 
gas. The very-hot plasma presents a severe energetics prob- 
lem, however: were it a hydrogen plasma its sound speed 
would be ~1500 km/s - considerably in excess of the local 
escape velocity of ~ 900 km/s ( Muno et al.[[20d4 ) - and it 
would escape (Yamauchi et al. 1990) on a short timescale. 



A steady-state situation would require a power to sustain 
this outflow considerably in excess of 10*" erg/s. A second 
difficulty is that there is no widely-accepted mechanism to 
heat the plasma to more than a few keV. 

There is no universally-accepted resolution to these 
anomalies. One interesting suggestion is that the 8 keV emis- 
sion is due to a very hot helium plasma which would be grav- 
itationally bound ( Belmont et al.|2005 ) . Another suggestion 
is that the 'plasma' is illusory, the emission actually being 
attributable to unresolved point sources. Recent Chandra 
observations around / = 0.08, b —1.42 (taken to be typical 
of the so-called X-ray Ridge) support this sort of picture 
( [Revnivtsev et a"L][2009[ ). On the other hand, the situation 
within the inner ~150 pc - where the 6.7 keV Fe line emis- 
sion strongly peaks ( Yamauchi et al.]|1990 l - may be quite 
different to that pertaining elsewhere in the Galaxy ( [Dogiel[ 



et al. 2009a). A deep observation (Muno et al. 2004) of the 



inner 17 with Chandra could explain only ^ 40% of the 
X-ray flux as due to dim point sources. Moreover, recent re- 
sults obtained with the SUZAKU X-ray telescope continue 
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to clearly suggest (Koyama et al. 2007 Dogiel et al.||20l6 l 
the existence of a hot plasma covering at least the central 
20'; this issue, therefore, has remained unresolved. 



A2 The GC magnetic field 

The amplitude of the large-scale magnetic field threading 
the GC region has remained uncertain for about twenty-five 
years (see Ferriere|2009 for a recent review) . Three different 
scales of field amplitude have been suggested in the light of 
different analyses: 



• mG field: ~GHz radio observations of the region reveal 
the presence of so-called "non-thermal filaments" (NTFs; 
see |Bicknell fc Li|[2001[ for a review of these fascinating 
structures), long and thin magnetic flux tubes that run pre- 
dominantly perpendicular to the Galactic plane through the 
central ~ 200 pc (diameter) and which are illuminated by 
synchrotron emission from relativistic electrons. The rigid- 
ity of the NTFs requires that they are suffused by mag- 
netic fields of ~mG amplitude (Yusef-Zadeh et al. 19871. 



A countervailing external pressure is then required lest the 
NTFs dissipate over short timescales given their large, inter- 
nal magnetic pressures and, given observational constraints, 
this can only be supplied by an external magnetic field also 
of ~ mG strength. This implies a ~ mG field pervades the 
region (| Morris fc Yusef-Zadeh|1989l ). 

• 100 /iG field: Alternatively, the assumption of pressure 
equilibrium between the various phases of the Galactic Cen- 
tre interstellar medium (including turbulent molecular gas; 
the contested ( [Revnivtsev et ciL]|2009[ ) "very hot" plasma; 
and the magnetic field) suggests fields of ~100 over ~400 
pc (diameter) size scales ( Spergel fc Blitz|[l992 L 

• 10 /iG field: Given their detection of the DNS source 
alluded to above, (LaRosa et al. 20051 have invoked the 



"equipartition" condition (minimizing the total energy in 
magnetic fields -I- cosmic rays, incidentally predicting rough 
equipartition in energy density between these two: |Burbidge| 
1956 Beck fc Krause 2005 1 to arrive at a large-scale field 
(over the central ~800 pc in diameter) amplitude of ~10 
within which cosmic ray electrons are synchrotron-radiating 
to produce the observed, non-thermal radio emission. 



2010a I 



As remarked above, some of us (Crocker et al, 
were recently able to show - with reference to the higher- 
frequency extension of the DNS radio spectrum in concert 
with GeV 7-ray data - that magnetic field amplitudes lower 
than 50 /xG are excluded over the central ~ 800 pc (in di- 
ameter). This work also calls into question the assumption 
of equipartition which may fail in the GC as well as in other 
star-forming regions (Thompson ct al. 2006). Other pieces 
of evidence - as described below - support the notion that 
the real DNS field is ~ 100 /xG. (This still implies an energy 
density more than 100 times larger than that typical for the 
field in the Galactic disk.) 



A2.1 Magnetic field required for equipartition with light 
field 

A magnetic field of ~ 60 /xG would be required were equipar- 
tition with the volumetric-average Uisnw — 90 eV cm~^ over 
the HESS region to be satisfied. In the Galactic plane - and 



many other astrophysical environments (see, e.g., Thompson 
et al.||2006 l- it can be inferred that the magnetic field and 
interstellar radiation field globally satisfy C/s ~ 2 f/isRF. 
A close-to-constant ratio between the two quantities cer- 
tainly seems to be necessary to maintain the linearity of the 



FRC (Lisenfeld et al. 1996 1. There are dynamic, local mecha- 



nisms that have been put forward to explain why this scaling 
should apply. These rely on the notion that magnetic fields 
arise in turbulent dynamo action with the ultimate source of 
such turbulence being the local, massive stellar population 
which injects turbulence via stellar winds and SN explosions 



(Lisenfeld et al. 19961 or, more directly, through radiation 
pressure (in starburst environments: Thompson|2008 1. 



A2.2 Magnetic field required for hydrostatic equilibrium 

A magnetic field exceeding that required to establish hydro- 
static equihbrium, B^qm, should lead to the development of 
the Parker ( 1966 1 instability and thus B^qm might be consid- 



ered an upper limit within the scenario explored by |Thomp-| 
son et al.l (120061) andlLacki et al.l (120091). The field required 



for hydrostatic equilibrium can be estimated as 2yjT,gaa'^tot 
mG where the gas and total surfa ce mass densities are given 
m g cm -2 ( [Thompson et al.|2006| ). The average surface mass 
density in gas over the HESS and DNS regions can be esti- 



mated on the basis of the work of Ferriere et al. ( 2007 1 to be 



0.07 g cm ^ and 0.01 g cm ^ respectively. The total surface 
mass density - dominated by stars - we infer from the work 



of Launhardt et al. ( 2002 \ to be around 21 times the gas sur- 



face density for the HESS region and 36 times for the DNS 
region. Hydrostatic equilibrium therefore requires fields of 
~ 600 /xG and ~ 200 ^.G for the HESS and DNS regions. 
It is noteworthy, then, that our previous work points to a 
field approaching this strength in the DNS. Interestingly, a 
Parker-like instability is suggested by mm-wave observations 
( Fu kui et al.|2006] ) of molecular filaments of several hundred 
parsec length within ~ 1 kpc of the GC (these measurements 
also suggest a ~ 100 /iG field amplitude). These observations 
have independently been interpreted as suggesting a ~ 100 
/iG field ( |Fukui et al.||2006| [Crocker et al.|2010a| . 



A2.3 Magnetic field in GC molecular clouds 

Measurements of the topology and amplitude of the mag- 
netic field interior to the GC's giant molecular cloud pop- 
ulation have been obtained via sub-millimetre polarimetry 



( Novak et al. 2003 Chuss et al. 2005 1 and Zeeman map- 



ping of the H I line in absorption (Crutcher et al. 1996) 



These measurements show that the cloud-interior fields are 
drawn out along a direction following the ridge of molecular 
and thermal radio emission, predominantly parallel with the 
Galactic plane. Such a configuration might arise naturally 
from the shearing action of the tidal field at the GC on the 
orbital motion on the molecular gas. Field amplitudes, more- 
over, reach ~mG levels within, e.g., the massive molecular 
envelope surrounding the dense Sgr B2 star-forming cores 
(Lis & Goldsmith[[1989| ICrutcher et al.[[1996[) with lower 



limits on larger scales at the 
1997l|Chuss et al.|2005 l. 



.^150 /iG level (^Novak et al 
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A2.4 Magnetic field topology 

The determination of the GC magnetic field's global topol- 
ogy is an area of active research. The geometry of the NTFs 
implies not only that the field is strong, but also that it is 
regular and poloidal ( Morris|2007 Ferriere|2009 and refer- 
ences therein). On the other hand, the field in dense molecu- 
lar material as revealed by the techniques mentioned above 
is predominantly parallel to the plane, at right angles to 
the large-scale poloidal field structure apparently traced by 
the NTFs. Most recently, NIR polarimetry measurements 
( Nishiyama et al.|2010 l - in principle, sensitive to magnetic 
field structure in both dense and rarefied ISM phases - have 
apparently revealed a global magnetic field structure that 
is parallel to the plane for Galactic latitudes jfoj < 0.4° but 
bends to run perpendicular at higher latitudes. 



A3 Volume of the source regions 

For the purposes of our modelling, motivated by the TeV and 
RC observations, for the HESS region we assume a cylindri- 
cal volume (seen in projection) which implies a volume of 
9.7 X 10®^ cm~"^. We assume the diffuse non-thermal source 
(DNS) emission volume to be an elliptical spheroid with cir- 
cular cross-section in the Galactic plane (i.e., its greatest 
line-of-sight extent is equal to its width along the plane). 
This implies a total volume of ~ 3.0 x 10®^ cm^. 



0.08 Me/yr for the DNS and HESS regions, respectively, 
and consistent with the upper-limit on the SFR within the 



inner 400 pc quoted by |Yusef-Zadeh et al.| (|2009j) of 0.15 
MQ/year (obtained by taking the stellar conte nt of the re- 



gion and an assumed age of 10 years; also see Schanne et 



al. (2007l). (Note that young, massive Main-Sequence stars 
are the predominant contributors to the luminosity of both 
Launhardt et al.|2002 |. 



regions: see 



B1.2 Interstellar radiation field 

We obtain the interstellar radiation fields (ISRFs) required 
by our modelling in the following fashion: i) for calculations 
concerning the DNS region, we employ the full interstel- 
lar radiation field determined for the Galaxy's inner 500 pc 
( Porter et al.|2006[ ). This has an energy density of ~ 19 eV 
cm~'^. It is to be noted that - as a check on self-consistency 
- taking the volumetric average of the radiation energy den- 
sity due to a point source of luminosity equal to Li^i'^ we 
find 17 eV cm~^, tolerably close to the result obtained from 
the work of [Porter et al.| ( |2006[ ) (obtained for a somewhat 
different volume). 

For ISRF of the HESS region, again using the Ltir 
of the region given above, we find a volumetric-average 
t/jgj^p — cm~^ over the HESS field; we assume the 

shape of the HESS ISRF spectrum is the same as that de- 
termined by [Porter et al.| ( |2006| l for the inner 500 pc of the 
Galaxy. 



APPENDIX B: STAR FORMATION AND 
SUPERNOVAE IN THE GC 

For an initial estimate of the star-formation rate (SFR) in 
the GC, we turn to the work of Figer et al. ( 2004 1. From this 
we infer that the HESS field is responsible for ~ 2% of the 



Galaxy's massive star formation ( Figer|2008 nominates the 
region inside r = 500 as responsible for 10% of the Galaxy's 
star formation). Given the SFR of the entire Galaxy is ~ 
4 Mo /year ( |Diehl et al.][2006| ), the HESS field SFR is ~ 
0.08 Mo/year. 



Bl GC star formation rate 

Bl.l Star-formation rate from FIR emission 

A more accurate handle on the SFR is potentially given by 
its FIR emission. On the basis of the work of iKennicuttI 
|T998l ) we have 

SFR = Mo/yr (Bl) 

2.2xl0«erg/s ^'^ ^ ' 

where the total infra-red (TIR) emission is given by Ltir — 
1.75 Lfir; Lfir = I/so Mm[l + Sioo pm/(2.58 Sao ^m)]; and 
L(d = 8kpc)6o Mm ^ 2.5 x 10*^ e rg/s Sao m„/(10^ Jansk y). 

On the basis of IRAS data ( [Launhardt et al.[[2002[ ) the 
Ltir of the DNS source region is 2.8 x 10 erg/s, and that of 
the HESS field 1.6 x lO^^ erg/s[^, implying SFRs of 0.12 and 



^ We note that the observed value of Ltir = 4.2 X 10* Lq for 
the HESS region is very close to the expectation generated by 
the total amount on H2 in the system (8.7 X 10^ Mq as inferred 
from the data of |Ferriere et ar][2007| l and the rough empirical 



B2 Supernova rate 

A first approximation to the supernova rate is given by the 
consideration that the HESS field is responsible for ~ 2% 
of the Galaxy's massive star formation and that the Galaxy 
experiences ~2 SN per century ( [Diehl et al.|2006| so that we 
expect around 0.04 SN per century in the HESS field. This 
simple estimate neglects effects like a possibly biased IMF 
in the GC; the detailed and independent arguments we set 
out below, however, confirm a central value for the SN rate 
of 0.04 per century with an uncertainty of a factor of ~2. 
It is significant for our purposes (as we explain below) that 
the estimate we derive for the SFR from the FIR output of 
the region (which traces the current SFR) is consistent with 
other measures. 



B2. 1 Supernova rate from FIR emission 



We infer on the basis of the work of Thompson et al. ( 2007 1 
that - as both the supernova rate and the total IR (TIR) 
luminosity are proportional to the star formation rate - one 
may relate the supernova rate to the Ltir via 



-^=8.9x10-" /JiT^ 
century Lq 



(B3) 



correlation between a given amount of (strictly, dense) H2 and 
the IR emission of a system established by |Mangum et aL[ l |2008l l: 



log 



2.0 -I- 1.0 log 



M 



H2 



-0 / V a/q 

on the basis of H2CO observations of starburst galaxies. 



(B2) 
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where /3i7 ~ 1, is described by Thompson et aL (20071. With 
Ltir ^ 4.2 X 10* Lq for the HESS field, we again determine 
a supernova rate of 0.04/century. 



B2.2 Supernova rate from stellar composition 

The type la supernova rate in the GC region defined by 
r < 250 pc and height < 50 pc has recently been estimated 
by Schanne et al. ( 2007 1 employing the scaling relations that 



have been determined ( Scannapieco fc Bildsten|2005 1 to re- 
late the supernova rate to the star formation rate, M and 
total stellar mass of external galaxies, M : 



RsNla{t) 

century 



M,{t) 
1010 Mq 



+ B 



M,{t) 



1010 MoGyr- 



(B4) 



where A = 0.044 ± 0.015 and B = 2.6 ± 1.1 ( [Schanne et al 
|2007[ ). Adopting now the stellar mass of the Galactic nuclear 



bulge (r < 300 pc) determined by Launhardt et al.| ( [2002 1 
of 1.4 X 10^ Mq, we can calculate the (central) value for the 
expectation of the SNla rate in this region of the Galaxy 
to be 0.04/century (where we also assume a formation time 
of 10 Gyr; Figer et al. 20041. For an estimate of the core- 



collapse SN rate, we incorporate the fact ( jMannucci et al.] 
20051) that the SN la rate is 0.35 ± 0.08 of the core-collapse 



rate in young stellar populations. Setting the first term in 
Eq. |B4| to zero, we determine that the expected core-collapse 
supernova rate in the nuclear bulge is 0.11/century (again 
taking central values) implying an overall rate of 0.15/cen- 
tury for all supernovae in the r < 250 pc region. 

This result must now be scaled to the region under con- 
sideration (somewhat smaller than the nuclear bulge as de- 
fined above). On the basis of the data presented in table 7 of 
( Launhardt et al.|2002p , there are two ways to do this: i) we 
can scale the above rate by the ratio of the inferred bolomet- 
ric luminosity of the HESS field to that of the larger r < 250 
pc region: 4 x 10*Lq/2.5 x 10^ Lq = 0.2; alternatively we 
can scale by the ratios of the inferred stellar masses of the 
two regions: 8 x 10** M0/I.4 x 10^ Mq = 0.6. This implies a 
total supernova rate in the HESS field of 0.02-0.08/century. 

Similar considerations to the above allow us to calculate 
that the SN rate within the entire DNS region (a volume ~30 
times larger than the HESS region) is only 1.5 — 2.0 times 
larger than in the HESS region alone. The SN activity is thus 
highly centrally peaked and, for simplicity, in our analysis we 
ignore the energy and particles injected by SNRs occurring 
in the DNS region but outside the HESS region. 



B2.3 SN rate upper limit from GC pulsar population 
The rates found above are consistent with that found by 



Lazio & Cordes ( 2008 1 and Deneva et al. ( 2009 1 from their 



radio studies of pulsars and radio point source pulsar can- 
didates in the inner 1° of the Galaxy, viz. that the GC SN 
rate has a rough upper bound at 0.1 per century (for a con- 
tinuous - rather than bursting - star formation history, as 
we favour: see below). 



B2.4 SN rate lower limit from X-ray emission 

We can arrive at a lower limit to the SN rate by requiring 
that there be enough power to support (only) the radia- 



tive loses of the plasma(s) and if we adopt the figure of 1% 
of SN kinetic energy lost into plasma heating from |Kaneda| 
|et al. (I997f. With the total X -ra y luminosity over the cen- 
tral region ~ 4 x 10^'^ erg/s (Belmont fc Tagger||2006| ), this 
consideration defines a lower limit of 0.01/century. 



B2.5 Gas turbulence 

The turbulent motions of the molecular gas throughout 
the region also represent a reservoir of energy which can 
be estimated to dissipate at a rate |Bradford ct al. 2005|) 
Lturb — 0A2vf^^i,/ Ad where Ad is a typical size scale for 
the turbulent structures. Estimating Vturb ~ 25 km/s and 
Ad ~ 10 pc we find Lturb ~ 4 X 10^^ erg/s over the ~ 9 x 10 
solar masses of H2 found in the HESS region. 



B2.6 High-velocity compact clouds 

A large-scale CO J=3-2 survey of the CMZ with the At- 



acama Submillimeter-wave Telescope Experiment (Oka et 



al. 20101 detects an unusual population of high-velocity 



{Av > 50 km/s) compact clouds (HVCCs) in the GC. The 
kinetic power of the HVCCs launched by explosive events 
in the HESS field is 3.8 x 10^^ erg/s, where the expansion 
velocity of each HVCC supplies the appropriate timescale to 
render a power from that HVCCs kinetic energy. 



B2.7 Combined constraint; thermalization efficiency 

Considering the power required to sustain both gas turbu- 
lence and the HVCCs in the region, we can place a lower 
limit on the supernova rate in the region of ~ 0.02 per 
century assuming ~ 100 % of the energy released per su- 
pernova goes into these two. Turning these considerations 
around, again with the requirement that both turbulence 
and HVCCs be sustained and adopting the maximum super- 
nova rate found above of 0.08/century we arrive at a conser- 
vative estimate of the minimum overall thermalization effi- 
ciency for GC supernovae of 77"ho"m = 0.4 (where, in this in- 
stance, thermalization refers to all power lost into heating or 
moving the GC ISM rather than radiativel^p^ . This is some- 
what larger than the ~10% thermalization efficiency deter- 
mined for normal or quiescent starforming galaxies (see, e.g., 
[Strickland fc Hcck man 2009), but this parameter is a func- 
tion of local ISM conditions with, in particular, higher gas 
densities tending to imply larger radiative losses (and, hence, 
smaZ/er thermalization values). It is then somewhat surpris- 
ing, a priori, that, in the molecular-gas rich milieu of the 
GC, the thermalization is rather high. In fact, further con- 
sideration reveals the high thermalization is plausible: super- 
novae explode into the low-density cavity created by their 



massive progenitor star's wind and UV output ( Chevalier 
1999). Moreover, in a starburst-like environment, one must 
consider the effect of previous SN which will tend to have 



1" Below, for purposes of determining the maximum and mini- 
mum star-formation-driven super-wind speeds we conservatively 
adopt r?J5ic?m = ^-^ ^^'^ '?thc"rm = whcrc, in this particular 
instance, we are concerned with only that proportion of the me- 
chanical power injected into the ISM that accelerates the overall 
outflow. 
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rendered the local ISM highly porous - riven with regions of 
high temperature, low density plasma - through their own 



gas heating (Hcckman et al. 1990 1. In any case, a thermaliza- 
tion efficiency of >40% is in the range recently determined 



by Strickland & Heckman ( 2009 1 for the M82 starburst 



B2. 8 Final estimate of SN rate 

We have, then, an estimate for the total supernova rate in 
the HESS field of 0.02-0.08/century, translating to a total 
kinetic power due to supernovae of LsN = l-Stol X 10""^ 
erg/s (adopting the fiducial 10^" erg per supernova into rel- 
ativistic particles decided above). 



The mass lost off massive stars is approximately equal 
to the mass returned to the ISM in supernovae and these 
components sum to give a total mass return of 

AMsN + AAf„,„d = AMtot ^ 0.02 - 0.03 Mg/year . (B5) 

This estimate is close to the range 0.03 — 0.1 M0/year esti- 
mated by [Figer et al.| ( |2004[ ) to be lost into a thermally- 
driven wind (over a somewhat larger region). Also note 
that, on the basis of §B1| the mass injected by super- 
novae and stellar winds per star formation rate is (0.02 — 
0.03) M© /year/ (0. 08 M© /year) = ( 0.25-0.38) close to the 
range of (0.19 — 0.31) arrived at by Strickland & Heckman 



(2009J in their detailed modelling of the starburst system 
M82. 



B3 Supernova mass injection 



We have previously estimated ( Crocker fc Aharonian|[20T0 1 
the hot mass injected into a potential superwind - that in- 
cludes both supernova ejecta and the swept-up mass - on 



the basis of the study of Zirakashvili & Volk ( |2006J of the 
nearby starburst NGC 253. On the basis of our previous 
study and adopting the central value of the supernova rate, 
0.04/century, one determines a mass injection rate due to 
supernovae in the HESS field of AMsn ^ 0.01 Mo/year. 

B4 Winds off massive stars 

Though expected to make only a subdominant (~10%) con- 
tribution to the total mechanical power starbursts in gen- 
eral, stellar winds may supply up to half of the mass re- 



turned to the ISM (Leitherer et al.| 


1992 


Leitherer & Heck- 


man 1995 Strickland & Heckman] 


2009 


1. This expectation 


is thought to be borne-out in the GC ( 


Muno et al.,,2004). 



Chevalier ( 1992 I originally calculated the wind power due 



to the massive (dominantly Wolf-Rayet) stars in the central 
pc of the Galaxy to be ~ 5 x 10'^'^ erg/s for a mass loss rate 
from He I emission line stars of 3 x 10""* Mo/year and an 
assumed wind velocity of 700 km/s. Over the entire central 5 
pc, the total mass loss rate from the Central s tellar cluster is 
around 3 x 10"^ Mg/year ( [Fryer et al.|2007 l and that from 
the Arches and Quintuplet clusters have been determined 
to be 4 X 10"* Mq /yeax and 3 x 10~^ M0/year, respectively 
( Lang et al.||2005 l. Given that Paschen-a studies (Wang et 
al.|2009 1 indicate that about half of the GC's massive stars 
are located outside these three clusters, we expect that stel- 
lar winds launch a total AMwind — 0.01 Mo/year. Assuming 
a typical wind speed of 700 km/s, this implies a total kinetic 
power of 1.3 x 10^^, consistent with the fiducial 10% of SN 
power expected. 



B6 Other sources of power 

For completeness, we mention that other sources of power 
considered in the GC context include i) pulsars, in partic- 
ular milli-second pulsars ( Wang et al.|200"6 l which are very 
difficult to detect in the GC environment because of pulse 
broadening of ~100 s at ~ GHz in the turbulent plasma of 



the GC (Cordes & Lazio 1997 Lazio & Cordes 20081 but 



which are generically expected to contribute around 10% of 
the kinetic power produced by SN; ii) tidal disruption of 
stars which approach too close to the central SMBH, a pro- 
cess predicted to occur at a frequency lO"**- 10~^ yr~^ , each 
event releasing io^^~^^ erg (Strubbe & Quataert 20091, pos- 



sibly leading to important radiative consequences ( Dogiel et 
|al.||2010| ); in) magnetic reconnection ( Tanuma et al.||1999[ ); 
iv) the gravitational energy of the region's giant molecular 
cloud tapped via the heating induced by the MHD waves 
excited by these clouds' passage through the (hypothesised) 



very hot (and highly viscous) plasma (Belmont & Tagger 
[2006| ). 



APPENDIX C: INDICATIONS AND 
EXPECTED CHARACTERISTICS OF A GC 
OUTFLOW 

There is considerable evidence - reviewed here - for a rather 
fast GC outfiow that would be expected to advect non- 
thermal particles from the system. Most of our discussion 
below will focus on a global GC outfiow beyond the central 
region (< 10 pc) surrounding Sgr A*. However, it should 
be noted that modelling of this environment (e.g.. Chevalier 
|1992[ [Ruffert fc Mefia||1994| | has itself provided compelling 
evidence for a substantial mass outfiow from the inner few 
lightyears of the Galaxy. 



B5 Total mass and power budgets 

Given the above, supernovae and (sub-dominantly) stellar 
winds inject, together, a mechanical power into the HESS 
region of ~ 1.4 x 10'"' erg/s, uncertain by a factor of ~2p^. 



There are a number of other sources of mechanical power in 
the system which, while ranging from the virtually-assured to the 
speculative, are all rather poorly constrained empirically. One of 
the results of our global analysis of the system's particle astro- 
physics is that there is no evidence that power above and be- 



Cl Direct evidence for a GC outflow 

The strongest evidence in support of the existence of an 
outfiow from the GC comes from observations of a ~ 1° (or 



yond that injected by supernovae (and stellar winds) is required 
in explaining the region's non-thermal emission at radio and TeV 
energies (the GeV spectrum may indeed, however, be dominated 
by the contribution of milli-second pulsars as discussed below). 
For completeness we list other potential sources of power below 
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~ 140 pc) tall and diameter < 130 pc shell of emission rising 
north of the Galactic plane called the Galactic Centre Lobe. 
This structure wEis first identified in 10 GHz radio continuum 
emission ( Sofue fc Handa|1984 1 and has been interpreted as 



evidence for a previous episode of either starburst ( Chevalier 



1992 Bland-Hawthorn & Cohen 2003 1 or nuclear activity 



(Kundt 1987 Melia fc Falcke||2001 ' or as due to magneto 



dynamic effects (Uchida et al. 11985 1. RC emission from the 



lobe's eastern part has HI absorption that clearly puts it in 
the GC region ( Lasenby et al.|1989 l. 

The coherent nature of the GC lobe was suggested 
by the discovery ( Bland-Hawthorn fc Cohen||2003 l at mid- 
infrared (MIR) wavelengths of filamentary structures coin- 
cident with the lobe's entire radio structur Modelling 
( Bland-Hawthorn fc Cohen|2003| ) of the MIR emission from 
shock-heated dust pointed to outflow masses and kinetic en- 
ergies of 5 X 10^ Mq and > 10^^ ergs, respectively, with an 
associated formation timescale of ~ 10^ years. 



Most recently Law et al. (20081 and Law (20101 have 



demonstrated the coherent nature of the lobe on the basis of 
its very similar morphologies in radio continuum, radio and 
optical recombination lines, and the MIR and |Law et al.| 
( 2009 1 demonstrated that its ionized gas is of high metallic- 



ity, consistent with a GC origin. This multi-wavelength anal- 
ysis shows, in addition, that the entire structure has a similar 
spectral index at ~GHz frequencies and that the structure 
is composed of multiple nested shells showing from, out- 
side to inside, respectively, i) MIR emission from entrained 
dust and poly-cyclic aromatic hydrocarbons, ii) synchrotron 
emission at RC wavelengths (diameter ~ 110 pc and with 
an equipartition (Beck & Krause"2005) field in the range 
40 - 100 /iG), and iii) optical and recombination line emis- 
sion from ionized material (which data also demonstrates 
the structure is located at the GC rather than being lo- 



cal). Also of note is that the detailed analysis of Law ( 2010 1 
considerably revises downwards (with respect to the esti- 
mate of Bland-Hawthorn fc Cohen||2003| the energetics of 
the structure to ~ 5 x lO'^'^ erg and finds that the formation 
of the Lobe is consistent with currently-observed pressures 
and rates of star-formation in the central few x 10 pc of the 
Galax>Pl 



C2 Hints from external galaxies 

In the case of external star-burst systems and active galaxies 



Veilleux et al. 



the evidence for outfiows is compelling (e.g. 
|2005[ ). Even for less luminous, star-forming systems akin the 
Milky Way (e.g., NGC 891) the evidence for outfiows from 



the central regions is strong ( Dahlem et al.|1994 Breitschw- 



lerdt 



2008[ ). For star formation surface densities exc eeding 



4 X 10 Mq yr kpc , empirical studies (Heck 



man et aL]|1990 Strickland fc Heckman|[2009 ) suggest that 



Another suggestive association is between the Radio Arc 
the brightest and most dramatic of the GC's non-thermal radio 
filaments — and the eastern extension of the lobe: 



(19841 



Yusef-Zadeh & Morris 



Pohl et al. 



Sotue k. Handal 



1 19921. 

lower limit on 



^(19881 

Note that even takmg the previously-claimec 
the total energy associated with the outflow of > 10^* erg, the 
mechanical power, 1.4 X 10-*° erg/s, supplied by GC supernovae 
occurring at the rate we have claimed can potentially supply this 
over the claimed formation timescale ~ lO'^ years. 



supernovae and massive stars act in concert to drive a mass- 
loaded 'super-wind' out of the starburst nucleus. (Alterna- 
tively, though not independently given the existence of the 
Schmidt star-formation law, galaxies that exceed a critical 



surface density of gas of 



0.05 g cm are observed to 



drive super- winds: Heckman et al.|2000 Lacki et al.||2009" l 

In general super- winds are comprised of multiple phases 
of varying density and velocity dStrickland fc Heckman 



|2009|). As revealed by observations of external, star-forming 
galaxies like M82, the 'wind fiuid' itself is a low-density and 
extremely hot plasma with an initial temperature in the 
range lO'' - 10*K. 

In any case, the areal star formation rate within the 
HESS region is ~2Mq yr'^ kpc-^ exceeding E^'' by 
almost two orders of magnitude. From the point of view 
furnished by studies of other galaxies, then, the expectation 
is that the GC should be driving an outflow. 



C2.1 Expected characteristics of a GC superwind 

As in external, star-forming galaxies, the mechanical power 
and mass injected collectively into the ISM by GC super- 
novae and massive stellar winds might be expected to gen- 
erate a 'super' wind with the following characteristics: 

(i) The asymptotic speed of such a wind (which the out- 
Volk et al.|[l990) should scale as 



flow quickly reaches: e.g. 

Vwind ~ \/ E/M which evaluates to 800-1000 km/s M range 
nominated above (0.02-0.03 Mq /year and a total mechani- 
cal power 011.4x10^" erg/sand assuming thermalization of 
100%). More accurately, the asymptotic speed of an outfiow 
powered by the thermalized energy injected by supernovae 
and other sources in the disk of a starburst galaxy can be 
rendered as (e.g., Zirakashvili fc Volk|2006 l : 



1200 km s~^^%^„ 



M 



E 

1.4 X 10*0 erg s-i 

1/2 



1/2 



0.03 Mq yr- 



(Cl) 



where E and M are the total energy and mass production 
rates and r^therm parameterizes the thermalization efficiency 
of the overall superwind. 

(ii) The central temperature of the wind-launching region 



can be rendered as ( Strickland fc Heckman||2009 1 



Te = 0.4 



E 



kB /3ic 
2 X lO'^ K 



M 



ria-i^ 



" 1/2 
E 



(C2) 



1.4 X 10« erg s" 



M 



0.03 Mq yr- 



where ^ is the mean particle mass (in units of the proton 
mass) normalized to the expected value for a H plasma and 
Aoad ^ 1 parameterizes the degree of mass-loading of the 
outflow. 

(iii) The central pressure of the wind-launching region can 
be rendered as ( Strickland fc Heckman|2009 1: 



0.118 



1/2 ol/2 



kB m 
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3.3 X 10'' K cm-^ vULpLUC 

^ 1/2 

X 



1.4 X 10*0 erg s-i 



M 

0.03 Mq yr- 



1/2 



65 pc 



(C3) 



where i?* is the radius of the starburst region which we have 
normahzed to the ~65 pc radius of the GC Lobe ( Law|2010 1 
(and consistent with the extent of the bright X-ray plasma 
region) and the ~ 1 is the participation factor defined 



by [Strickland fc Heckman ( 2009 1 which parameterizes the 



fraction of total star formation that is actually involved in 
driving the superwind. 

(iv) Lastly the central density of the wind-launching re- 
gion can be rendered as ( [Strickland &: Heckman|2009[ ): 



= 0.296 



,1/2 



i;i/2 Ri 



E 



-1/2 ^3/2 ^ 
-1/2 



1.4 X 10« erg s-i 



M 



0.03 Mq yr- 



3/2 



65 pc 



(C4) 



C3 



Additional evidence in support of the 
existence of GC outflow 



C3.1 Gamma-ray deficit in inner Galaxy 

A long-standing problem in CR studies ( [Strong fc Mattox| 
1996[ [Breitschwerdt et al.[[2002 1 concerns why the 7-ray in- 
tensity across the inner Galaxy increases (relative to local) 
so much more slowly (for decreasing Galactic radius) than 
expected given the increase in areal density of plausible CR 
sources (and also given the growth in the density of ambient 
gas) towards the centre of the Galaxy. 



This has a natural explanation (Suchkov et al. 1993 



Breitschwerdt et al. 2002 1 if particle escape times also de- 



cline in proportion to the increase in source density and this, 
in turn, is naturally explained if the same objects that are 
the ultimate sources of power for the non-thermal particles 
also serve to drive an increasingly fast wind to advect par- 



ticles out of the Galactic plane (also cf. Blitz et al. 1993 1 



Following this logic, GC CR density should roughly satisfy: 



ncR oc — — (C5) 

With ~2% of the Galaxy's star formation occurring in the 
inner 300 pc (diameter), the over-abundance of CRs inside 
300 pc with respect to the value in the Galactic disk should 
then be 



''CR 



0.02 



15 kpc 
150 pc 



20 



wind 

f "wind 
I 1000 km/s 



where the (diffusive) escape of CRs fro m the disk is parame- 
terized via an effective v'^'"'' = 10 km/s ( [Suchkov et al.|l993[ ) 
and we normalize to a wind speed of order the gravitational 



escape speed (as observed for external galaxies: e.g. 
et al.|2009| |. 



Heesen 



Implicitly here we find additional evidence that, in fact, 
a wind of a least a few 100 km/s is required in the GC: 
7-ray studies at ~GeV energies of the GC that failed to 
detect emission associated with some of the massive molec- 
ular clouds (Sagittarius B, C, or D) belonging to the CMZ 
( |Mayer-Hasselwander et al.[ [1998') require that the GC CR 
energy density be no more than a few times larger than the 
local value in the few GeV energy range. Note, however, 
that this constraint requires that CRs in this energy range 
sample all the mass of the ambient molecular clouds, includ- 
ing, in particular, the material in the cores of the clouds - 
and theoretical pogiel fc Sharov[[1990[ [Gabici et al.|[2007 l 



and e mpirical studies ([Protheroe et al. 
2009[F1 



2008 Fujita et al. 



show this may not be true, potentially weakening 



this argument. 



C3.2 Evidence of GC outflow on larger scales 

Across the electromagnetic spectrum, there is considerable 
evidence of GC outflows on scales much larger than the ~ 1° 
extension of the GC lobe, up to order tens of degrees. This 
evidence includes: 

(i) The North Polar Spur (NFS), a thermal X-ray/radio 
loop extending ~ 80° north of the Galactic plane first iden- 
tified by [Sofue[ ( [T977t ; [Sofue fc Handa[ ( [T984| ; [Sofue[ ( [2000| ) 
and interpreted by them as evidence for a ^ lO'^" erg nu- 
clear outburst ~ 15 Myr ago. Note, however, that there are 
also strong arguments to suggest the NFS is a local ISM fea- 
ture produced by the continuous effects of supernovae and 
stellar winds from the ScoCen association (see, e.g., [Mill er 
et al.|2008| |. 

(ii) On scales of up to ~ 25° north of the GC, the Galactic 



centre spur identified ( Sofue et al. 1989 1 in 408 MHz ( Haslam 
[et aHl982) and 1408 MHz ( [Reich et al.|1990p data. 

(iii) In general, ROSAT X-ray observations covering the 



central regions (inner few kpcs) of the Galajcy ( Sofue"2000| 
[ Almy et al 12000] [Yao fc Wang|2007||Everett et al. 2008; Bre^ 
'itschwerdt "2008^ seem to indicate a nuclear or bulge outflow 
of plasma with a temperature few x 10® K. The detailed 
analysis of [Everett et al.[ ( [2008[ ) suggested a Galactic wind 
towards the centre of the Galaxy asymptoting to 760 km/s in 
the halo and strongly suggested a role for CRs in helping to 
launch this wind (interestingly, however, these authors find 
that - assuming constant ISM conditions fixed to match the 
local ones - a wind cannot be launched inside ~ 1.5 kpc due 
to the growth in the Galaxy's gravitational potential). 

(iv) Radio polarimetry measurements (Duncan et al 
1998| which show a plume of emission extending ~ 15° 
north of the Galactic plane, roughly centred on the GC, 
with weaker evidence for a mirroring southern structure. 

(v) From 



Lockman ( 1984 1 deter- 



21 cm observations, 
mined that the inner Galaxy seems to be relatively depleted 
in H I, plausibly attributed to clearing of the neutral gas 
from the region \z\ > 500 pc, R ^ 3 kpc by a nuclear out- 
flow. 



Also see Jones, Crocker et al. 2010, 'Australia Telescope Com- 
pact Array Radio Continuum 1384 and 2368 MHz Observations 
of Sgr B', submitted to the Astronomical Journal. 
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(vi) The infrared dust observations by Bland-Hawthorn 



& Cohen ( 2003 \ indicate an outflow structure on very large 
scales, approaching the size of the North Polar Spur. 

(vii) A number of studies find evidence (in UV absorp- 
tion features of various background objects) of high- velocity 
material plausibly associated with a large-scale, Galactic nu- 



clear outflow. Keeney et al. ( 2006 1 find evidence (in sight- 



lines towards two high-latitude active galactic nuclei) of ma- 
terial both escaping from and falling toward the centre of the 
Galaxy with a typical speed of 250 km/s, the latter feasibly 
indicating the fall-back of material in a bound Galactic wind 
or 'fountain' (see below). Zech et al. ( 2008[) report the detec- 



tion of high velocity clouds (HVC) of, significantly, super- 
solar metallicity along the sightline of star at a distance ~ 
7.5 kpc. In their view, the probable close proximity of the 
HVC material to the GC, its complex, multiphase structure, 
its high metallicity, and its large velocity all indicate an as- 
sociation with a Galactic nuclear outflow or fountain. 

(viii) Though the connection to the GC is more circum- 
stantial, it has long been recognized that a flow of hot gas 
from the Galactic disk to the halo is necessary to sustain the 
Shull|1996 l. 



latter (see, e.g., 



Prime candidates for the requi- 
site gas conduits are the so-called chimneys of gas that have 
been observed to form above the 'super-bubbles' around OB 
associations in external galaxies and also, apparently, the 
Galactic disk ( Normandeau et al. [1996 1. Given the massive 
stellar population of the GC mirrors that found in such asso- 
ciations, it seems not unlikely that, by analogy, a chimney- 
like structure should form above the GC especially given the 
detection of super-bubble like structures associated with the 



GC massive stellar clusters (see Rodriguez-Fernandez et al, 
2001 and references therein). 



Finally, we note the existence of other less-direct indi- 
cations of a GC outflow include the following: 



i) One explanation ( Shore & LaRosa 1999 1 of the non- 



thermal radio filaments (NTFs) posits that they are, in fact, 
akin to cometary plasma tails formed from the interaction of 
a hypothesised, large-scale, magnetized, ~1000 km/s Galac- 
tic wind 'draping' dense molecular clouds in the GC region. 



(ii) Lu et al. ( 2008 1 report the detection of 10 X-ray fila- 
ments in the inner ~ 15 pc around GC which they interpret 
as pulsar wind nebulae whose tails, mostly pointing out- 
wards, again suggest a GC wind of velocity ^ 400 km/s. 



Markoff (2010) summarizes evidence from Chandra for bi- 
polar lobes of hot gas extending over similar scales around 
Sgr A*. Mid-infrared observations of two cometary struc- 
tures on arcsecond scales around Sgr A* ( Muzic et al.|2007 1 
have recently been interpreted ( |Muzic et al.||2010[ ) as iiidi- 
cating the presence of a coUimated outflow of at least a few 
100 km/s driven either by the SMBH or stars in its vicinity. 

(iii) Extreme Doppler-broadening of the 1.8 MeV line 
from "^^ Al decay is evident in emission detected around the 
GC shows and indicates the decaying nuclei have a (non- 
thermal) velocity of ~500 km/s ( |Naya et al.||1996| , consis- 
tent with a wind advecting the material out of the Galactic 
plane. 

C4 Stalling of GC outflow 

For an adiabatic outflow powered by a wind, a sufficient con- 
dition for the flow to reach inflnity is that the asymptotic 
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Table Dl. Observational data used to derive the spectrum for 
the HESS and DNS regions identified in the main text. Notes: 
*At 74 MHz the large-angle Galactic plane synchrotron back- 
ground/foreground flux contribution is not measured (and hence 
not accounted for) due to the interferometric nature of the 
VLA observation. tTotal flux (measured with Greenbank single 
dish) supplied by Dr Crystal Brogan (private communication). 
** Greenbank Earth Station fluxes are lower limits. Key: VLA 
= Very Large Array, GBT = Greenbank telescope, EBG = Ef- 
fclsbcrg, PKS = Parkes, GBT(ES) = Greenbank Earth Station, 
NOB = Nobeyama. 



wind velocity be faster than the gravitational escape speed. 
Given that the Galactic escape speed on relevant scales is 
~900 km/s, this condition will be satisfied over at least some 
of the parameter space we favour. However, for a radiative 
fiow, the conditions on the wind velocity are more stringent 
if it is not to stall. The case of a wind driven from a super 
stellar cluster located close to the GC (e.g., the Arches clus- 



ter) has recently been explored by Rodriguez-Gonzalez et al 
(2009). These authors determine the stall height of an out- 
fiow as a function of metallicity and wind velocity and find 
over the relevant parameter space (500 km/s< v^ind < 800 
and 1 < Z < 10) that heights between 1 and 50 kpc from the 
plane can be achieved (note that radio recombination line 
observations have found an unusually low electron temper- 
ature that implies a high metal abundance for the ionized 
gas in the GC lobe: [Law et aL||2009 l. 



APPENDIX D: RADIO DATA 

Sources of radio data for the HESS and DNS radio spectra 
are listed in table iDll 

In addition to the radio data employed in the paper, in 
table [Dl we also list 8.35 and 14.35 GHz data (Langston 
|2000| ). Unfortunately, the data at these frequencies 



et al 



were processed with a median filtering algorithm that arti- 
ficially attenuates the fiux density of sources larger than ~ 
34' so that these fiux densities strictly represent lower lim- 
its. Given the strong concentration of the 8.35 GHz emission 
into structures on size scales less than ~ 0.5°, however, we 
treat the 8.35 GHz fiux as a datum to be fit by our proce- 



dure. The 14.35 MHz datum shows the expected ( Finkbeiner 
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et al. 2004 1 up-break above 10 GHz due to the rising con- 



tribution from spinning dust and is, therefore, not adopted 
as a constraint in our fitting. 



APPENDIX E: DETAILED JUSTIFICATION OF 
ASSUMPTION THAT PARTICLE TRANSPORT 
IS ENERGY-INDEPENDENT 

In our modelhng we assume that the non-thermal particle 
transport timescale is energy-independent, a situation that 
is realised if, for instance, the removal of particles from the 
region is affected by their being advected away on a wind 
(rather than diffusing in an energy-dependent fashion by 
scattering on magnetic turbulence). This implies that trans- 
port does not alter the slope of the injection spectrum (of 
either relativistic electrons or protons over the energy range 
where its associated timescale is the least). Also note that - 
as the pp cross-section is, to a fair approximation, also energy 
independent - the steady-state, in situ proton spectrum has 
the same shape as the injection distribution (cf. the situa- 
tion for the steady-state, in situ electron spectrum where the 
different energy dependences of the various electron cooling 
processes certainly can introduce spectral features). 

Although an energy-independent timescale is an as- 
sumption, we believe this to be physically well-motivated 
for a number of reasons 

(i) There is no spectral evidence for diffusion processing 
of the spectrum. The spectral index of the diffuse ~ TeV 
emission is 2.2-2.3. It is firstly difficult to reconcile such 
a hard spectrum with any diffusion processing having oc- 
curred; in general, the least steepening, A7 = — 1/3 occurs 
for a Kolmogorov spectrum of turbulence. In this case, the 
injection spectral index must be very hard, 1.9-2.0. (Any 
other spectrum of turbulence would require an even harder 
injection spectrum.) Given that - as we explain at length 
below - we expect that supernova remnants accelerate the 
protons and primary electrons needed to explain the non- 
thermal emission detected from the region and that the lo- 
cal CR spectrum (also due mostly to acceleration by SNRs 
in the Galactic disk) has a spectral index ^ 2.7, the most 
natural explanation is that both GC and local CRs have the 
same spectrum at injection (in the same class of astrophys- 
ical object) but, whereas the local CRs have been diffusion 
steepened, the GC CRs have not. The only escape clauses 
to this argument would seem to be: 

(a) the system is out of steady-state with insufficient 
time having elapsed since an assumed, single injection 
event for the particle distribution to have been steepened 



(cf. the interpretation adopted by Aharonian et al.|2006 
that a single CR-injection event ~ 10 year ago at the GC 
explains the observed diminution in the 7-ray to molecular 
column ratio beyond \l\ ~ 1°). As argued above, however, 
all the evidence is consistent with the system being in or 
close to steady-state. 



(b) the diffusive transport time, tdffsn 



7(4D), 



actually close to the rectilinear escape time trect ~ r/c 
for a region size r. In this case, the diffusive escape time 
must become energy-independent so that the steady-state, 
in situ spectrum must flatten back to that at injection 
( Aharonian fc Neronov|2005 1. An exploration of this sce- 



nario is beyond the scope of the present work but note, 
at least as far as the protons are concerned, this situation 
still realises an effectively energy-independent transport 
time as we have modelled. For an application of this idea 
to the GeV/TeV point source coincident with Sgr A* see 
Chernyakova et al. (20101. 



A related point is that, within the admittedly large statisti- 
cal errors, the spectral index of the diffuse ~ TeV emission 
is constant over the solid angle of the emission region. There 
is thus no evidence for steepening of the CR hadron distri- 
bution as the particles are transported - and age - away 
from their putative discrete sources. 

(ii) The radio spectral information is also highly sugges- 
tive of energy-independent transport in the HESS field. Ac- 
cording to the modelling results of |Lisenfeld fc Volk| ( |2000| ), 
radio spectral indices a ^ 0.7 (for 5*1^ oc h'~") indicate such 
and a power-law fit to the radio data from the HESS field 
for all radio data 1.4-10 GHz finds a spectral index of 
a ^ 0.54, corresponding to an emitting electron population 
(for Fe oc E~'^) with spectral index 7 ~ 2(a -I- 1) ^ 2.1 (in- 
clusion of the 74 MHz datum leads to an even harder spectra 
but this has been affected by free-free absorption). 

(iii) Diffusive transport has been shown to be very slow in 
the GC by direct modelling of particle trajectories ( |Wommer| 
et al.|2008| ). 

(iv) All the evidence from starbursts - to which we have 
argued the GC is closely analogous - is that CR transport 



is overwhelmingly advective (Suchkov et al. 1993 Dogiel 



& Breitschwerdt 20091. Theoretically, this is expected in 



light of the fact that CR scattering produced by self-excited 
Alfven waves should be very effective in starbursts, leading 
to marked reductions in the predicted diffusion coefficient 
( [Ptuskin et al.|1997| [Zirakashvih fc Volk|2006| ). On the em- 
pirical side, in the well-known star-burst system M82, for 
instance, cosmic ray escape is probably dominated by ad- 
vection out of the gas midplane in a starburst-driven wind 
(jThomps on et al.|[2007l [Klein et alT]|1988| [Seaquist fc Ode- 
|gard„199lf 



